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Abstract 

Following Symanzik we argue that the Schrodinger functional in lattice 
gauge theories without matter fields has a well-defined continuum limit. Due 
to gauge invariance no extra counter terms are required. The Schrodinger 
functional is, moreover, accessible to numerical simulations. It may hence 
be used to study the scaling properties of the theory and in particular the 
evolution of the renormalized gauge coupling from low to high energies. A 
concrete proposition along this line is made and the necessary perturbative 
analysis of the Schrodinger functional is carried through to 1-loop order. 



1. Introduction 



The Schrodinger functional is the propagation kernel for going from some 
field configuration at time x'^ = to some other configuration at = T. In eu- 
clidean space-time it can be written as a functional integral over all fields with 
the specified initial and final values. The renormalization of the Schrodinger 
functional has been discussed by Symanzik [1] some time ago in the course of 
his proof of the existence of the Schrodinger picture in renormalizable quan- 
tum field theories (for an introduction to this paper see ref.[2]). His result is 
that the Schrodinger functional can be renormalized by adding the usual coun- 
terterms to the action plus possibly a set of further terms that are integrals of 
local polynomials in the field and its derivatives over the x° = and a;° = T 
hyper-planes. In the case of the 0'* theory, for example, two polynomials, 0^ 
and (l>dQ(j), are needed. Pure gauge theories are simpler in this respect, because 
no extra counterterms are required here, as we shall argue in sect. 2. If we 
choose a lattice to regularize the theory, this statement simply means that the 
Schrodinger functional converges in the continuum limit, provided, of course, 
the bare coupling is scaled in the usual way. 

Our motivation to consider the Schrodinger functional is that we would 
like to apply the finite size scaling technique of ref.[3] to gauge theories. The 
ultimate goal of this programme is to compute the running coupling in say 
the minimal subtraction (MS) scheme of dimensional regularization at short 
distances given in units of the low energy scales of the theory (the string 
tension, for example, or the mass of the lightest glueball). In other words, 
our aim is to connect the non-perturbative infrared behaviour of the theory 
with the high energy regime, where the coupling is logarithmically decreasing 
according to the renormalization group. 

To explain in which way the Schrodinger functional enters this calculation, 
we need to recall the basic strategy of ref.[3]. One begins by putting the 
theory in a finite spatial volume with linear extension L and periodic boundary 
conditions in all directions. Next one introduces some renormalized coupling 
(f'{L) which docs not depend on any scale other than L and which can hence 
be considered a running coupling. This coupling is then computed over a 
range of L through numerical simulation of the lattice theory, using a recursive 
procedure which allows one to go from large values of L (where contact is made 
with the non-perturbative scales) to the perturbativc small L domain. At 
these distances g^{L) can be analytically related to other more commonly used 
couplings such as the coupling in the MS scheme of dimensional regularization. 
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Note that the latter are usually defined in infinite volume. All reference to a 
finite volume thus disappears from the final result. 

The precise definition of the running coupling g^{L) is not of principal 
importance. The coupling should however be accurately computable through 
numerical simulation and its scaling properties should not be strongly influ- 
enced by the presence of a non-zero lattice spacing. It is our experience that 
these requirements are not easy to fulfil. In particular, extracting a running 
coupling from correlation functions of Wilson loops is difficult, because large 
loops give a poor signal while small loops are affected by lattice artefacts. As 
we shall explain in sect. 2, running couplings can be straightforwardly defined 
through the Schrodinger functional. With carefully chosen boundary values 
for the gauge field, these couplings have, moreover, the desired technical prop- 
erties which make our finite size scaling study feasible. 

The idea to probe gauge theories through external fields is not new of 
course. In particular, the background field method [4-8] has long proved to 
be an efficient tool for perturbative calculations, in a wide range of theories 
and contexts (see e.g. refs.[9-12]). Our work can also be understood as a 
continuation of rcf.[13]. Most of the basic concepts that wc now exploit already 
appear there, in a language oriented towards statistical mechanics. 

In the present paper the Schrodinger functional in Yang-Mills theories is 
discussed in detail. Numerical results on the scaling behaviour of the couplings 
derived from it will be published elsewhere. After going through the formal 
definition of the functional in sect. 2, we verify explicitly, using dimensional 
regularization, that indeed there are no additional counterterms at one-loop 
order of perturbation theory (sect. 3). We then discuss the definition and 
cutoff dependence of the Schrodinger functional on the lattice (sect. 4). In our 
numerical work, we take constant Abelian fields for the boundary values of 
the gauge field at x° = and x° = T, and so, in sects. 5 and 6, we set up the 
perturbation expansion for this case. In particular, the stability of the gauge 
field configuration around which one expands must be established and the 
gauge fixing must be done with care to obtain all boundary terms correctly. 
As an application we calculate our running coupling to one-loop order, for the 
theory with gauge group SU(2) (sect. 7). 
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2. A first look at the Schrodinger functional 



The aim in this section is to give an introduction to the Schrodinger 
functional in Yang-Mills theories and to establish the basic notation. We 
shall not at this point worry about the mathematical status of the quantities 
considered. Later on when the theory is regularized either dimensionally or by 
passing to the lattice formulation, we shall be able to deal with the Schrodinger 
functional on a more rigorous level. 

For definiteness wc assume that the gauge group is SU(A^). Our conven- 
tions regarding indices, group generators etc. are listed in appendix A. 

2. 1 Formal definition of the Schrodinger functional 

Our starting point is the Hamiltonian formulation of the theory in the temporal 
gauge. In this framework the theory is specified at a fixed time, say = 0, 
by assuming canonical commutation relations among the basic field variables 
and by giving the Hamilton operator. 

As explained in the introduction, the Schrodinger functional will be used 
to study the scaling properties of the theory in finite volume, and so we take 
space to be an L X L X L box with periodic boundary conditions. SU(A/^) gauge 
fields are accordingly represented by periodic vector potentials ^/c(x) on IR"^ 
with values in the Lie algebra of SU(A'') f . To preserve periodicity under gauge 
transformations 

Afc(x) ^ At(x) = A(x)Afc(x)A(x)-i + A(x)5feA(x)-\ (2.1) 

only periodic gauge functions A(x) will be admitted. A can thus be regarded 
as a mapping from a 3-dimensional torus to S'\J{N). Continuous functions of 
this kind are topologically non-trivial in general. More precisely, they fall into 
disconnected topological classes labelled by an integer winding number 

^ = ^£ d^x ekij tr { {AdkA-') (AS^A'^) (AO,- A'l) } . (2.2) 

Gauge transformations with non-zero winding number are explicitly allowed 
and will play an important role later on. 

f We do not consider the possibility of twisted periodic boundary conditions in this 
paper. 
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In the Schrodinger representation the quantum mechanical states of the 
theory arc wave functionals ip[A], where A runs over all gauge fields as de- 
scribed above. A scalar product is formally given by 

(V^IX) = / B[A] i^iAYxiA], B[A] = dA^(x). (2.3) 

x,fc,a 

Only gauge invariant states il^[A], i.e. those satisfying 

^P[A'^]=^[A] (2.4) 

for all gauge transformations A, are physical. In particular, we choose the 
vacuum angle 6 [15,16] to vanish. Any given wave functional tp[A] can be 
projected on the physical subspace through 

i;[A] ^ Pi;[A] = j D[A] i;[A% D[A] = ndA(x), (2.5) 

where dC/, U S SU(A^), denotes the normalized invariant measure on SU(A^). 

The gauge field A'^ (x) can in the obvious way be interpreted as an operator 
field acting on wave functionals ip[A]. The canonically conjugate field is the 
colour electric field 

The magnetic components of the colour field tensor are 

F,",(x) = dkAti^) - diAli^) + r'"=4(x)^?(x), (2.7) 
and the Hamilton operator IH is given by 

IH = j^d'x ||F«,(x)Fo"fe(x) + ^F,«,(x)F,«,(x)| (2.8) 

with 5o being the (bare) gauge coupling. 

For any smooth classical gauge field Cfc(x), a state \C) may be introduced 
such that 

(CIV') = ^l^[C] (2.9) 



for all wave functionals ijj[A]. \C) is not gauge invariant, of course, but it 
may be made so by applying the projector IP. The (euclidean) Schrodinger 
functional Z[C', C] is now defined by 

Z[C',C] = {C'\e-'^^\P\C). (2.10) 

We shall always assume T > and do not explicitly indicate the dependence 
of the Schrodinger functional on this parameter. If we insert an orthonormal 
basis IV'n), n = 0, 1, 2, . . ., of gauge invariant energy eigenstates, the spectral 
representation 

oo 

Z[C', C] = J2 e-^-^'^nlC'mcr (2.11) 

n=0 

is obtained, where En are the energy eigenvalues (the spectrum is discrete 
in finite volume). Since only physical intermediate states contribute, it is 
evident that Z[C',C] is invariant under arbitrary gauge transformations of 
the boundary fields C and C 

2.2 Functional integral representation 

The matrix elements of the euclidean time evolution operator e~"^"^ between 
gauge invariant states can be expressed through a functional integral over 
all gauge field configurations Afj_{x) in four dimensions with < x° < T and 
periodic boundary conditions in the spatial directions. For the matrix element 
(2.10), the appropriate initial and final values of the gauge field are 

r C^(x) at x° = 0, 
AJx) = (2.12) 

and an integration over all gauge transformations A is required to account for 
the projector IP [cf. eq.(2.5)]. The functional integral representation of the 
Schrodinger functional thus reads 

Z[C', C]= J D[A] J B[A] e-^[^], (2.13) 

where the measure D[j4] now stands for an integration over all components of 
the euclidean field. The euclidean action is given by 

S[A] = /" d^xtr {F^.F^.} , (2.14) 
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with 

i^M- = d^A^ - d^A^ + [A^, A^]. (2.15) 

Finally, an overall normalization factor independent of the boundary values C 

and C has been dropped in eq.(2.13). In the following we do not keep track of 
such factors, because we are ultimately only interested in the field dependence 
of the Schrodinger functional. 

The reader may wonder at this point how it comes that the time compo- 
nent of the gauge field appears in the functional integral, while in the Hamil- 
tonian formulation it did not. The reason is that the integral (2.13) and the 
boundary conditions (2.12) are invariant under the gauge transformation 

A^{x) n{x)A^{x)n{x)-'^ + n{x)d^,n{x)-'^, (2.16) 

A(x)^Q(x)Uo=oA(x), (2.17) 

provided the gauge function Q,(x) satisfies 

n{x) = 1 at x° = T. (2.18) 

An admissible gauge fixing condition for this symmetry is the temporal gauge 
Aq = 0. The associated Faddeev-Popov determinant is field independent and 
may be included in the overall normalization of the integral. The time compo- 
nent of the gauge field can thus be eliminated and it is now straightforward to 
make contact with the Hamiltonian expression for the Schrodinger fTinctional. 

For the perturbative calculations to be discussed later on, other gauge 
choices will be more convenient so that in the following we keep the time 
component of the gauge field as an integration variable. 

We may however use the gauge invariance of the functional integral to 
reduce the integration over the gauge transformation A to a sum over topolog- 
ical classes. To this end, first note that after the inner integral in eq.(2.13) has 
been performed, one is left with some function of A to be integrated over. This 
function depends only on the winding number of A, because any two gauge 
functions in the same topological class may be connected through a gauge 
transformation Q as described above [cqs.(2.16-18)]. So if we choose, for each 
integer n, some fixed representative A„(x) in the class of gauge functions with 
winding number n, the Schrodinger functional becomes 

00 „ 

Z[C',C]= ^ /D[^]e-^[^], (2.19) 

n=— 00 
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where now we require that 

f C^"(x) at = 0, 
Ak(x) = \ ^ ^ ' (2.20) 

It is convenient to set Aq = 1, but an exphcit choice of the non-trivial A„'s 
would not be useful in what follows. That a sum over topological classes is 
needed on top of the functional integration is perhaps not too surprising, if 
we recall that this is the normal situation in Yang-Mills theories on compact 
manifolds without boundary. 

2.3 Instanton bound 

It is well-known [14] that the gauge field action is bounded by 

S[A]>^\Q[A]\, (2.21) 
9o 

where 

Q[A] = -r^ f d^xiv{F^,*F^,} (2.22) 
denotes the topological charge of A and 

Pp,u — '2^IJ,i'pcr^p(y (2.23) 

the dual of the field tensor (2.15). Q[A\ may be expressed through the bound- 
ary values C and C and the winding number n. To see this note that 

|tr {F^/F^.} = e^^p^d^iv {A,dpA„ + lA^ApA^] . (2.24) 

Taking the boundary conditions (2.20) into account and integrating by parts 
then yields (after some algebra) 

Q[A\ = Scs[C']-Scs[C\-n, (2.25) 

where 

Scs[C\ = J d^xeuj tr {C\diCj + ICkQCj} (2.26) 
is the Chern-Simons action of the boundary field C. 
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2-4 Induced background field 

For small couplings the functional integral (2.19) is dominated by the fields 
around the absolute minima of the action. From the instanton bound it follows 

that we only need to inspect a small number of topological sectors to find the 
fields with least action (and the specified boundary values). In general there 
are several gauge inequivalent minimal action configurations and these can 
even occur in different winding number sectors. 

For simplicity wc shall from now on restrict attention to boundary values 
C and C , where the absolute minimum is attained in the n = sector and 
where, furthermore, the minimal action configuration B^{x) is unique up to 
gauge transformations. This is the typical situation if C and C" are small. In 
the following B will be referred to as the induced background field. 

For given boundary values C and C , it is usually impossible to obtain the 
induced background field in closed analytical form. Wc may, however, invert 
the procedure and start from any known solution B of the field equations. If 
we define C and C through 

Cfc(x) = Sfc(x)|,o=o ' C'ki-^) = Bk{x)Lo=T , (2-27) 

the boundary conditions will be trivially satisfied. We are then left with 
the problem to prove that the chosen field B is the unique minimal action 
configuration with these boundary values. 

Little work is needed to show this, if the field tensor 

G^. = d^B, - d^B, + [B^, B,] (2.28) 

is self-dual, 

G/xi, = *Gnvi (2.29) 

and if, furthermore, the bound 

Scs[C']-Scs[0\<\ (2.30) 

holds. The point is that the instanton bound (2.21) is saturated in this case. 
All fields A in the n = sector hence have an action greater or equal to 

S[B] = ^ {ScslC] - Scs[C]} . (2.31) 
9o 
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And the same is true for all other sectors, as one may quickly show by com- 
bining cqs. (2. 21), (2.25), (2.31) with the bound (2.30). 

Wc have thus proved that under the conditions stated above, B is an 
absolute minimum of the action. It is in fact the only one, up to gauge 
transformations, because any other minimizing configuration would have to 
be self-dual and satisfy the same boundary conditions. Since the self-duality 
equation is first order in the time derivatives, the solution in the temporal 
gauge is uniquely determined by the initial values at x*^ = 0. In a general 
gauge, this means that any two solutions are related by a gauge transformation. 

A simple example of a self-dual field is obtained by making the ansatz 

Boix) = 0, Bk{x) = b{x^)Ik. (2.32) 

Here, 6 is a real function and the group generators Ik are chosen such that 
they form an irreducible representation of the angular momentum algebra 

[Ik,li]=ekijlj. (2.33) 

Up to unitary transformations, there exists only one irreducible representation 
of this algebra with dimension N. In particular, the square of the angular 
momentum is given by 

44 = -\{N^ - !)• (2.34) 

Note that the action of the field (2.32) is finite if 6(x°) is smooth and if L < oo. 
With the ansatz (2.32) the self-duality equation (2.29) reduces to 

dob = b'^, (2.35) 

so that 

6(a;0) = (r - x°)-^ (2.36) 

To ensure the regularity of the solution in the interval < x*^ < T, the 
integration constant r must either be negative or greater than T. We leave it 
to the reader to work out the condition (2.30), but it is obviously satisfied for 
sufficiently large values of r. Wc have thus found a one-parameter family of 
globally stable background fields. 
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2. 5 Renormalization of the Schrddinger functional 

In the weak coupling domain, the Schrodinger functional can be computed by 
performing a saddle point expansion of the functional integral (2.19) about 
the induced background field B. For the effective action 

T[B] = -In Z[C',C], (2.37) 

an asymptotic series of the form 

T[B]=g^^ro[B] + Ti[B]+g^o^2[B] + ... (2.38) 

is thus obtained f. The leading term is given by 

To[B] = giS[B], (2.39) 

while the higher order contributions are sums of vacuum bubble Feynnian dia- 
grams with an increasing number of loops. The Feynman rules involve vertices 
and propagators that depend on the background field. To make the diagrams 
well-defined an ultra-violet regularization will be needed at this point. We 
shall either use dimensional regularization or introduce a space-time lattice. 
In both cases the gauge invariance of the theory is preserved, which is crucial 
for the renormalization to work out in the way described below. 

Initially the saddle point expansion of the Schrodinger functional is per- 
formed in the regularized bare theory. That is, we integrate over the bare 
gauge field A, use the bare action to generate the Feynman rules and impose 
the boundary conditions (2.20) on the bare field. Since the regularization 
respects the gauge symmetry, the effective action r[B] does not depend on 
the gauge fixing condition employed. Moreover, it is invariant under gauge 
transformations 

B^{x) B^{x) = n{x)B^{x)n{x)-^ + n{x)d^n{x)-^ (2.40) 

of the background field for arbitrary (periodic) gauge functions ft. 

Ultra-violet divergences now appear in each order of the bare coupling go 
when we try to remove the regularization. We certainly expect that some of 

f As already mentioned earlier, we are only interested in the field dependence of the 
Schrodinger functional. Any additive contributions to the effective action independent of B 
are hence dropped without further notice. 
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the divergent terms are cancelled by the usual renormalization of the coupling 
constant. As we shall show in detail in the following sections, the one- loop 
effective action actually becomes finite after the coupling is renormalized. In 
other words, what is suggested is that the (bare) Schrodinger functional does 
not need any renormalization besides the coupling constant renormalization. 

A rigorous proof of this conjecture to all orders of perturbation theory 
is beyond the scope of this paper. We may, however, make our proposition 
more plausible by recalling Symanzik's work on the Schrodinger representa- 
tion in quantum field theory [1]. Symanzik studied the massless ^'^ theory 
and showed that all divergences in the Schrodinger functional can be can- 
celled by renormalizing the coupling constant and by including the boundary 
counterterms 

/ d^x {Z^cf)^ + Z24>do(p} + [ d^x - Z2(/)5o0} (2.41) 

in the bare action. The terms proportional to (fP' do not influence the prop- 
agators and vertices and just add to the effective action. Because the renor- 
malization constant Zi is linearly divergent, they are not needed if one em- 
ploys dimensional rcgularization. The other counterterms are equivalent to 
a rescaling of the boundary values of the field [1]. They are logarithmically 
divergent and must be taken into account when deriving the renormalization 
group equation for the Schrodinger functional. 

Symanzik also expressed the expectation that in a general renormalizable 
field theory, the Schrodinger functional can be made finite by renormalizing the 
coupling and mass parameters and by adding a few boundary counterterms to 
the action. These are proportional to some local composite fields of dimension 
less than or equal to 3 integrated over the .x" = and =T hypcr-plancs. 

Coming back to Yang- Mills theories, the obvious question is whether there 
are any candidates for such counterterms. Since the effective action is indepen- 
dent of the gauge fixing condition employed, we should be able to write them 
in a gauge invariant form not involving the Faddeev-Popov ghosts. We now 
note, however, that any non-trivial gauge invariant polynomial in the gauge 
potential and its derivatives has dimension greater than 3. A countcrterm pro- 
portional to the Chern-Simons action (2.26) is also excluded, because the latter 
is odd under parity. The upshot then is that divergent boundary terms cannot 
occur in this theory and that consequently the (bare) Schrodinger functional 
is finite after the coupling has been renormalized. 
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2. 6 Running couplings in finite volume 

As explained in sect. 1, the finite size scaling study that we propose to apply 
to Yang-Mills theories is based on the idea of a renormalized coupling g^ (L) , 

which runs with the box size L. Starting from the Schrodinger functional, 
there are many ways to introduce such a coupling. For example, wc may 
choose some background field B which depends on a dimensionless parameter 
rj. Prom the above we then infer that 

r'[s] = (2.42) 

is a renormalization group invariant. A renormalized coupling may hence be 
defined through 

f=r'o[B]/r'[B] (2.43) 

[cf. eq.(2.38)]. Note that Tq[B] is just a normalizing factor, which guarantees 
that g^ coincides with Qq to leading order of perturbation theory. 

In general the chosen background field depends on several external scales 
and the corresponding renormalized coupling must be regarded as a function 
of all these parameters. In addition it depends on the box size L and the 
propagation time T. To obtain a coupling which runs with L, we simply scale 
all dimensionful parameters with fixed proportions relative to L. For example, 
we may set T = L and take the self-dual configuration (2. 32), (2. 36) with 

r = -L/ri (2.44) 

as the background field. 

There are, of course, other choices of background fields and a correspond- 
ing manifold of renormalized couplings. In particular, a simple alternative are 
spatially constant Abelian fields, which, for reasons given later, are actually 
more suitable for our numerical work. 



12 



3. Perturbation Expansion 



We now discuss the saddle point expansion of the Schrodinger functional 
around the induced background field using dimensional regularization. In 
particular, we would like to show that to one-loop order the rcnormalization 
of the Schrodinger functional works out in the way described above. Our 
calculations rely on well-established techniques, which have previously been 
applied in the context of instantons and the semi-classical approximation [17- 
22]. There is also a formal similarity with the so-called background field 
method [4-8], although here we do not expand in powers of the background 
field. 

As before we shall assume that the absolute minimum of the action occurs 
in the n = sector and that the minimal action configuration B is unique up to 
gauge transformations. To avoid some technicalities when fixing the gauge, we 
shall in addition require that the gauge group acts freely on the boundary val- 
ues C and C. This means that the only spatial gauge transformations A which 
leave C or C invariant are constant and equal to a central element of SU(A''). 
The self-dual background field considered earlier has this property. Abelian 
fields however do not and the gauge fixing procedure must consequently be 
reexamined in this case (see sect. 6). 

3.1 Dimensional regularization 

The familiar techniques of dimensional regularization apply to Feynman dia- 
grams in the momentum space representation [23-26]. A different approach 
is required in the presence of background fields, because the propagators and 
vertices do not in general have their standard form. One possibility is to stay 
in position space and to insert the heat kernel representation for the propa- 
gators. The Feynman integrals then become well-defined analytic functions of 
the space-time dimension. In the following we shall proceed along these lines 
and refer the reader to ref.[27] for an introduction to the method. 

Dimensional regularization starts by extending space-time to a D dimen- 
sional manifold, where the extra D — 4 dimensions are here assumed to be 
spatial with the usual periodic boundary conditions. The bare action of a 
gauge field An{x) on this manifold is again given by eq.(2.14), except that 
we now integrate over all D coordinates of x and that the Lorentz indices 
run from to D — 1 (cf. appendix A). The boundary values C and C and 
the background field B are taken to be independent of the extra coordinates. 
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Their components in these directions are set to zero. In particular, 



S[B]=L^-'{S[B]}^^„ (3.1) 

and the same volume factor also appears in all higher order contributions to 
the effective action, because the translation symmetry in the added dimensions 
is not affected by the background field. 

The renormalization of the gauge coupling is performed as usual. We 
choose the minimal subtraction (MS) scheme [28] and accordingly denote the 
renormalized coupling by 5ms- The relation between the bare and renormal- 
ized couplings reads 

5o = m''5ms { ^ + £ (^)^^s I , D = A-2e, (3.2) 

where ^ is the normalization mass. The singular coefficients 

I 

fc=i 

coincide with the counterterms calculated in ordinary perturbation theory. In 
particular, the one-loop coefficient is given by 

and the two- and three-loop coefficients can be found in refs. [29-31]. 
3.2 Gauge fixing 

The action S[A] and the a priori measure D[^] are invariant under arbitrary 
gauge transformations A [cf. eq.(2.40)]. We are interested in evaluating 

the functional integral for fixed boundary values C and C , and so, in this 
subsection, restrict attention to the group Q of transformations $7 which leave 
the boundary values intact. Our assumptions on the boundary values imply 
that a gauge function $7 belongs to Q if and only if 

f e'^'^™/^ at x° = 0, 

I gi27rm /N ^0 ^ ^ 
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for some integers m and m'. Q accordingly decomposes into iV^ disconnected 
components. 

The center of Q consists of all gauge functions Q = e«27rm/Ar^ These 
transformations act trivially on gauge fields and so it is really the factor group 
Q = Q/Zn with which we are concerned when fixing the gauge. Q may be 
identified with the m' = component of ^. It acts freely on the space of 
gauge fields, i.e. = A for some A implies = 1. 

The gauge fixing procedure is the usual one: we add a gauge fixing term 
to the action and include an integration over Faddeev-Popov ghost fields c and 
c with the appropriate action. To write down the gauge fixing term, the gauge 
field A integrated over is decomposed according to 

A^{x) = B^{x) + goq^ix). (3.6) 

The "quantum" field q introduced here is the new integration variable, while 
B will be kept fixed. In perturbation theory, only the winding number zero 
sector contributes. The boundary conditions (2.20) thus become 

qk{x) = at x° = and = T. (3.7) 

Note that the time component qo remains unrestricted at this point. 

In terms of the quantum field q, a gauge transformation A — > A^ amounts 
to the substitution 

q^ix) ^ n{x)q^{x)n{x)-' + g^^ [B^{x) - B^(x)] . (3.8) 

In particular, if we set 

n{x) = l-goio{x) + 0{gl), (3.9) 
the transformation becomes 

q^{x)^q^{x)+D^uj{x) + 0{go), D^ = d^+AdB^. (3.10) 
This suggests that we take 

S^i[B,q] = -Ao / d'^xtr {D^q^D^q,} (3.11) 
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as the gauge fixing term, where Aq denotes the bare gauge fixing parameter. 
This gauge is referred to as the background gauge. The associated Faddeev- 
Popov action reads 



Spp[B,q,c,c] = 2 J d"" X tr {c D ^{D^+ go Ad q^)c}, (3.12) 
and so we end up with 

e-r[B] ^ J Qj^j J D[c]D[c]e-^*°-'[^'«''='^], (3.13) 



where 



StotAB, q, c, c] = S[B + goq] + Sgi[B, q] + Sfp[B, q, c, c] (3.14) 

is the total action. 

We have aheady noted above that the spatial components of q must vanish 
at x*^ = and = T. The boundary conditions on qo and the Faddeev-Popov 
fields c and c are determined by the gauge fixing procedure. In particular, they 
depend on the choice of the gauge fixing term and on the gauge transformation 
properties of the boundary values C and C. The influence of the latter is 
subtle and will only become clear after a while. 

In the background gauge and for boundary values as specified at the 
beginning of this section, it turns out that 

Doqo{x) = c{x) = c{x) = at x° = and = T (3.15) 

are the correct boundary conditions. This is not easy to prove, because the 
functional integral is not really well-defined before we fix the gauge. A rigorous 
derivation can (and will) be given after we pass to the lattice formulation of 
the theory (cf. sect. 6). 

At this point it is still possible to explain, on a heuristic level, how the 
boundary conditions (3.15) arise. Let us first consider the Faddeev-Popov 
fields. When going through the gauge fixing procedure, one notes that the 
Faddeev-Popov operator 

Afp = -D^ {Df,+goAdq^) (3.16) 

acts in the linear space of infinitesimal gauge transformations. So if t is an 
anti-commuting parameter, the transformation n(x) = exp[tc{x)] must be an 
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element of Q, and the same is true for c. Prom eq.(3.5) we now infer that the 
ghost fields vanish at x'^ = and = T. 

The boundary conditions on qq are harder to justify. An obvious case is 
the background field Landau gauge D^qn = 0, where the vanishing of Doqo 
at the boundaries is an immediate consequence of the boundary conditions 
on Qf.. For general values of the gauge fixing parameter, we may resort to 
the slightly obscure argument that the gauge fixing function Df^qj^ should be 
a mapping from the space of gauge fields to the space of infinitesimal gauge 
transformations. Since the latter vanish at x*^ = and x*^ = T, we again 
conclude that qo satisfies the boundary conditions (3.15). 

3.3 One-loop effective action 

It is now straightforward to expand the total action and the gauge fixed func- 
tional integral (3.13) in powers of go. To second order we have 

StotAB,q,c,c] = S[B]-2 J d^xtr {ig^ Ak?^ + c Aqc} + 0(50), (3.17) 

where the operators Ai and Aq are defined through 

Aig^ = -D,D,q^ + (1 - \o)D^D,q, - 2[G^,,g,], (3.18) 
Aqc = -D^D^c. (3.19) 

Ai and Aq act on fields in D dimensions and are to be distinguished from 
Ai = A 



and Aq = Aq 

D=4 



(3.20) 

D=4 



Both are hermitean operators relative to the obvious scalar products in the 
spaces of fields satisfying the boundary conditions (3.7) and (3.15). They are 
also elliptic (if Aq > 0) and bounded from below. In particular, there exists a 
complete set of eigenfunctions with a discrete spectrum of real eigenvalues. 

From the definition (3.19) it follows that all eigenvalues of Aq are strictly 
positive. Negative eigenvalues of Ai are excluded too, because 5 is a minimal 
action configuration. There could be a finite number of zero modes, but this 
is an unlikely case which we shall not discuss any further, i.e. our assumptions 
on B are supplemented accordingly. 

If we now insert eq.(3.17) in the functional integral and expand the inte- 
grand, we obtain the series (2.38) for the effective action. The leading term 
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is given by eq.(2.39) as expected. At one-loop order we only need to perform 
the Gaussian integrals over the quantum field and the Faddeev-Popov ghosts 
and so end up with 

ri[S] = I Indet Ai -Indet Aq. (3.21) 

The determinants occuring here are defined according to the rules of dimen- 
sional regularization (cf. sect. 3.3 of ref.[27]). 

In the following it is convenient to set Aq = 1- We are free to make 
this choice, because Ti[B] is independent of the gauge fixing parameter. A 
rigorous proof of this statement could actually be given at this point, but since 
the argument is a bit lengthy (and the result is expected anyway), we do not 
present it here. 

A more explicit expression for the determinants, exhibiting the depen- 
dence on the space-time dimensionality, can be given in terms of the heat 
kernels of Ai and Aq. Heat kernels have simple factorization properties, and 
if we take these into account, one gets [27] 

IndetAi = -J J (iVe-*^ ) (TVe"*^^ - 2£TVe-*^°) , (3.22) 

lndetAo = -y ^ (Xre-*^ ) Tre"*^". (3.23) 

The first factor in these integrals arises from the added dimensions. It involves 
the Laplacian A^ on the circle and is explicitly given by 

1 °° 

TVe-*^' = J2 e-*(2Wi')\ (3.24) 

n= — 00 

The other factors, Trc"*^^ and Trc~*^°, are defined in four dimensions and 
are thus independent of e. Finally, the term proportional to e in eq.(3.22) must 
be included because the quantum field q has D vector components rather than 
4. 

An important point to note is that the integrals over the proper time 
t are convergent for Ree > 2. At the upper end of the integration range, 
convergence is guaranteed by the exponential fall-off of the integrand (Ai and 
Aq are positive). For t ^ 0, on the other hand, the asymptotic behaviour 
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of the heat kernels is computable. This goes under the name of the Seeley- 
DeWitt expansion and will be discussed below. The result is that the functions 
integrated over are proportional to t~^~^^ at small t and so are integrable for 
Re£ > 2. 

3.4 Seeley-DeWitt expansion 

To evaluate the one-loop effective action near four dimensions, the integrals 
(3.22) and (3.23) must be analytically continued from Ree > 2 to a region 
including the origin e = 0. The key result which enables us to do so is the 
Seeley-DeWitt expansion. There exists an extensive literature on the subject 
(an incomplete list of references is [32-41]) and we shall, therefore, be rather 
brief. 

Let A stand for one of the operators Ai or Aq. The Seeley-DeWitt 
expansion states that 

Tre-*^ "2^"^^ + "3/2(A) t'^/^ + ai(A) r ^ + . . . , (3.25) 

where the coefficients aj/2{/^) are algebraically computable. There are two 
different kinds of contributions, the volume and the boundary terms. Vol- 
ume terms only occur for even j and arc equal to a local polynomial in the 
background field B and its derivatives integrated over the space-time mani- 
fold. The boundary terms are constructed similarly, except that the integral 
is taken over the hyper-planes at a;° = and x'^ = T (with possibly a sign 
difference between the two). The dimension of the polynomials must be equal 
to 4 — J and 3 — j, respectively. 

A crucial observation now is that the spectrum of Ai is invariant under 
the substitution B — > B^, for any gauge function $7. We simply transform the 
eigenf unctions covariantly, viz. 

q^{x) ^ n{x)q^{x)n{x)-\ (3.26) 

and then note that the eigenvalue equation and the boundary conditions are 
preserved. The same remark also applies to Aq and so it is clear that the 
coefficients aj/2{^) must be gauge invariant. 

The gauge symmetry just described is not accidental. The complete gauge 
fixed functional integral is actually invariant if B is gauge transformed and the 
fields integrated over are rotated covariantly. In particular, as we have noted 
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earlier, the effective action satisfies 

r[B^] = r[B]. (3.27) 

We emphasize that this property does not depend on our choice of gauge fixing 
term. The effective action is independent of the latter. The advantage of the 
background gauge only is that every individual contribution to the effective 
action is already invariant, while in a general gauge this would not be the case. 

Besides the trivial polynomial (the constant) there is no gauge invariant 
polynomial in the background field and its derivatives with dimension less than 
4. The polynomials with dimension 4 are linear combinations of tr {Gjui^Gpo-}. 
Taking this into account it is clear that all coefficients aj/2iA) with j > are 
independent of the background field. Furthermore, using one of the known 
computational techniques (ref.[39] for example), one finds that 

AT [■ 

ao(Ao) = -^ao(Ai) = ^ j d^xtx {G^.G^.} . (3.28) 

It has certainly not escaped the reader's notice that our argumentation here 
parallels the discussion of the renormalization of the Schrodinger functional in 
subsect. 2.5. In particular, it is because of gauge invariance that we are able 
to exclude the presence of boundary terms in the coefficients ao(A). As we 
shall see below, this is intimately related to the fact that the renormalization 
of the effective action does not require any extra counterterms. 

3.5 Removal of the ultra-violet cutoff 

To perform the analytic continuation of the integrals (3.22) and (3.23) towards 
e = 0, it is useful to break up the integration range at t = 1. Since the 
integrands fall off exponentially, it is clear that the integrals from 1 to infinity 
are entire analytic functions of e. So we only need to care about the integration 
from to 1. In that range the integrands may be split into singular and regular 
parts according to 

4 

TVe"*^ = ^aj/2(A)t-^/2^i?(t|A). (3.29) 

From the Seeley-DeWitt expansion we know that i?(t|A) is of order t^^"^ at 
small t. Furthermore, using the Poisson summation formula, one may show 
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that 



Tre 



L 



(47rt)i/2 ^ 

^ ^ r) = — rv 



(3.30) 



i.e. this factor is proportional to The contributions of the regular parts 

to the integrals are hence analytic in the region Ree > — |. The remaining 
terms are easy to work out, since the integrands are explicitly known. 

To be able to write the result of these computations concisely, we intro- 
duce the zeta function 

C(s|A) = TrA-^ (3.31) 

The trace converges for Re s > 2, but after passing to the heat kernel repre- 
sentation, one may show, following the steps taken above, that C(*|A) extends 
to a meromorphic function in the whole complex plane. In particular, there is 
no singularity at s = and we may define 



c'(o|A) = ^c(^ia: 



(3.32) 



s=0 



The determinants of Ai and Aq are now given by 



In det Ai 
In det Aq 



- + In 47r - 7e + ]^ 



L-2^ao(Ai)-C'(0|Ai), (3.33) 



— h In 47r — 7e 



L-2-ao(Ao)-C'(0|Ao), 



(3.34) 



where 7e = 0.577... denotes Euler's constant and all terms vanishing at e = 
have been dropped. 

Prom these results we immediately deduce that 

This singularity is exactly cancelled by the coupling constant renormalization. 
Indeed, from eq.(3.2) we get 

V[B] = / J- - 2i(£)| Vo[B] + T,[B] + 0(5L), (3-36) 
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and after inserting eq.(3.4), we are left with 



iC'(0|Ai) + C'(0|Ao) + 0(5^s) 



(3.37) 



We have thus shown that to one-loop order the Schrodinger functional renor- 
malizes in the expected way. 

The zeta functions appearing in eq.(3.37) are complicated functionals of 
the background field, which cannot normally be reduced to simpler expres- 
sions. For spatially constant background fields they can be worked out nu- 
merically, essentially by listing all eigenvalues of Ai and Aq up to a certain 
size (cf. sect. 7). 



4. Lattice Formulation 



A non-perturbatively meaningful definition of the Schrodinger functional 
can be given in the framework of lattice gauge theories. The lattice regu- 
larization is not unique — a choice of lattice action must be made and some 
further arbitrariness is involved when the Schrodinger functional is introduced. 
We expect that these details do not matter in the continuum limit and what 
follows should therefore be regarded as one possible approach to the problem. 

4-1 Gauge fields 

We choose to set up the lattice theory in euclidean space and thus imagine that 
a regular hypercubic lattice is superimposed on the space-time manifold. T is 
assumed to be an integer multiple of the lattice spacing a so that the possible 
values of the time coordinate of a lattice point x arc 0, a, 2a, . . . ,T. 

A lattice gauge field U is an assignment of a link variable U {x, ji) € SU(A^) 
to every pair (x, x + afi) of lattice points (/t denotes the unit vector in the ji- 
direction). In particular, the temporal link variables U{x,^) are defined for 
all lattice points with < x° < T. Gauge functions 0(a;) live on the sites x of 
the lattice and take values in SU(A/"). They act on the link variables according 
to 

U{x, ii) U^{x, ii) = n{x)U{x, ii)n{x + aft)-^. (4.1) 
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As in the continuum we require that gauge fields and gauge transformation 
functions are periodic in spatial directions with period L (which must also be 
an integer multiple of a). 

The lattice regularized Schrodinger functional will be defined later on as 
an integral over all lattice gauge fields with prescribed boundary values 

W{x,k)= U{x,k)\^o=o and W\K,k) = U{x,k)lo=T ■ M 

To make contact with the continuum definition of the Schrodinger functional, 
W and W must be related to the continuum boundary values C and C 

A natural relationship is suggested if wc recall that U{x, n) is the parallel 
transporter for colour vectors from x + afi to x along the straight line con- 
necting these two points. To achieve the best possible matching of lattice and 
continuum boundary values, we should thus identify VF(x, k) with the corre- 
sponding parallel transporter determined by the field Cfe(x). In other words, 
we set ^ 

W(y:,k) =Vexp!^a dt Cfc(x + ak - tak)| , (4.3) 

and W'{x,k) is similarly given by the field C^(x) [in eq.(4.3) the symbol V 
implies a path ordered exponential such that the fields at the larger values 
of the integration variable t come first]. Note that this construction is gauge 
covariant: if we perform a gauge transformation C C^, the associated 
boundary field W transforms as a lattice gauge field should. 

4-2 Lattice action and the Schrodinger functional 

Following Wilson [42] the action S[U] of a lattice gauge field U is taken to be 

^[^] = iE^(f)t'-{l-^(f)}' (4.4) 

% p 

where the sum runs over all oriented plaqucttcs p on the lattice and U{p) 
denotes the parallel transporter around p. The weight w{p) is equal to 1 in all 
cases except for the spatial plaquettes at x° = and x^ = T which are given 
the weight ^. The significance of this weight factor will be discussed later on. 
The lattice regularized Schrodinger functional is now defined through 

Z[C',C] = I D[i7]e-^[^1, D[U]=l[dU{x,,,), (4.5) 

X,IJ. 
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where one integrates over all lattice gauge fields U with fixed boundary values 
as described above [recall that dU denotes the normalized invariant measure 
on SU(A^); the product in eq.(4.5) is over all lattice points x and directions 
such that the link (x, x+ap,) is in the interior of the lattice] . We emphasize that 
the lattice Schrodinger functional is regarded as a functional of the continuous 
fields C and C rather than the boundary values of the link field. The latter 
are determined through eq.(4.3) and its primed relative. In the continuum 
limit, C and C are kept fixed, while the lattice spacing a is sent to zero and 
the bare coupling is scaled according to the renormalization group. 

Compared to the formal continuum expression (2.19) we seem to miss a 
sum over topological classes here. Such an average is not needed on the lat- 
tice, because the functional integral (4.5) is already invariant under arbitrary 
gauge transformations of the boundary fields. To sec this, recall that a gauge 
transformation C — > induces a corresponding transformation W of 
the lattice boundary field at x° = 0. In the functional integral, such a change 
of the boundary values can be compensated by a substitution U of the 

integration variables, where is an arbitrary gauge transformation function 
with 

, , r A(x) at 2:° = 0, , , 

n{x) = \ ^ ^ ' (4.6) 

I 1 at = T. 

The crucial point to note is that gauge functions with these boundary values 
exist on the lattice, independently of whether A is topologically trivial or not. 

A quantum mechanical interpretation of the lattice Schrodinger functional 
(4.5) can be given through the well-known transfer matrix construction [43- 
47]. The transfer matrix Tq in the temporal gauge f/(x, 0) = 1 is a hermitean 
operator which acts on Schrodinger wave functions at time x° = 0. It 

may be regarded as the lattice expression for the step evolution operator e~"^" 
f. As in the continuum theory, physical wave functions iplW] must be gauge 
invariant. The associated projector is again denoted by IP and we may also 
introduce a state \W) in a way analogous to |C) [cf. eq.(2.9)]. 

The point we wish to make is that our definition of the lattice Schrodinger 
functional is precisely such that the identity 

Z[C', C] = {W'\ (To)^/" P\W) (4.7) 

f In some of the papers quoted above it is the gauge projected operator IT = TToP 
which is referred to as the transfer matrix. There is no difference between IT and Tq on the 
physical subspace. We here stick to Tq to match with the notation employed in sect. 2. 
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holds. This formula is entirely analogous to the corresponding continuum 
expression (2.10). In particular, a spectral representation of the type (2.11) 
may be obtained, for all values of the lattice spacing. 

4-3 Background field 

For given boundary values C and C (with properties as specified earlier) and 
sufficiently small lattice spacings, we expect that the absolute minimum of the 
lattice action S[U] is non-dcgcncratc up to gauge transformations. The mini- 
mizing configuration V should moreover converge to the induced background 
field B of the continuum theory, 

V{x,fi) = l + aB^{x) + Oia^), (4.8) 

provided that both fields are transformed to a definite gauge (the temporal 
gauge, for example). In perturbation theory, the lattice background field V 
plays the same role as B did in the continuum. It is, therefore, often necessary 
to find an analytical or at least an accurate numerical representation of V. 

To illustrate these remarks, let us consider the self-dual background field 
which was introduced in sect. 2 [eqs.(2.32)-(2.36)]. Since C and C" are spa- 
tially constant in this case, the lattice boundary fields are simply given by 

W{x, k) = exp {ab{0)lk} and W'{x, k) = exp {ab{T)Ik} . (4.9) 

The associated minimal action configuration V must be a solution of the field 
equations which one derives from the Wilson action S[U]. To be able to write 
them concisely, we introduce the plaquette field 

P{x, fi, u) = V{x, fi)V{x + afi, i')V{x + az>, ii)~^V{x, (4.10) 

and its covariant divergence 

d*P{x, ii) = 

3 

{P{x, fj,, v) — V{x — au, u)~^P{x — av, fx, i')V{x — au, i^)} . (4.11) 

i/=0 

One may now show that the lattice action is stationary if and only if the 
traceless anti-hermitean part of d*P{x,ii) vanishes, 

d*P{x,fi)-d*P{x,^)^ - ^tT{d*P{x,n)-d*P{x,n)''} = 0, (4.12) 
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for all links {x,x + afi) in the interior of the lattice. 

The form of the continuum background field suggests that we try to solve 
the lattice field equations (4.12) with the ansatz 

V{x,(}) = l, V{x,k) =ex-p{av{x°)Ik} , (4.13) 

where v is some real function to be determined. In the case of the SU(2) 
theory, a solution is in fact obtained in this way, if v satisfies 

3 

0<x° <T. (4.14) 

The operator d in this equation denotes the forward lattice derivative, 

a/(xO) = i[/(xO + a)-/(xO)], (4.15) 

and d* the backward derivative. The ansatz (4.13) also works out for SU(3), 
but for > 4 wc sTispcct that a more complicated expression is needed, with 
at least two unknown functions. 

In the practically relevant range of T/a, say T/a = 10, ... , 100, the so- 
lution of cq.(4.14) with the required boundary values (and the lowest action) 
can be determined numerically to any desired precision. Alternatively, if we 
assume that an expansion of the form 

i;(x°) = voix'^) + avi{x^) + 0^2(2;°) + . . . , (4.16) 

holds, where Vk{x^) are smooth functions independent of a, the resulting tower 
of equations, 

divo = 2vl (4.17) 
d^vi = Gv^vi, (4.18) 

etc., can be solved iteratively. We expect, of course, that the leading term vq 
coincides with the continuum solution b. This is consistent with the lowest 
order equation (4.17), and from eq.(4.18) we now infer that vi is a linear 
combination of and 6""^. If we take into account that vi must vanish at 
x° = and = T, a little calculation shows that vi is in fact equal to zero 
everywhere, and so we conclude that 

v{x°) = b{x°) + 0{a''). (4.19) 



d* 



— sin (|a^5v) = 8 cos {^av) 



— sin (hav) 

a ^■^ > 
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We have also worked out the order correction and compared our result with 
the numerical solution of eq.(4.14). Complete agreement was found, and the 
approach of the lattice solution to the continuum background field is, therefore, 
well under control. 

One point, however, that remains to be checked is that the lattice back- 
ground field so constructed is indeed a configuration with least action. To con- 
vince oneself that there are no other lower minima, one may run a relaxation 
program, on a range of lattices, starting from various initial configurations. In 
the case of most interest to us, the constant Abelian fields, such a numerical 
"proof" of stability is fortunately not needed, because the minimum property 
of these fields can be established by analytical means (cf. sect. 5). 

4.4 Continuum limit in perturbation theory 

The effective action on the lattice, T[B], is defined through eq.(2.37) and is 
considered to be a functional of the continuum background field. We are 
certainly free to do so, because the Schrodinger functional (4.5) depends on 
the boundary values C and C and these are in one-to-one correspondence 
with B. The notation is also in accord with our understanding of the lattice 
as a device to regularize the theory which is to be removed at fixed B. 

The perturbation expansion (2.38) of the lattice effective action may be 
derived following the steps taken in sect. 3 for the case of dimensional regular- 
ization. In particular, the leading term ro[i3] is equal to Qq times the action 
of the lattice background field V. An interesting new aspect is that the gauge 
fixing procedure can be carried out rigorously. This will be discussed in sect. 6 
and an explicit one-loop calculation is presented in sect. 7. Our aim here is to 
describe how the continuum limit is reached in perturbation theory. 

Symanzik's analysis [48,49] of the cutoff dependence of Feynman diagrams 
on the lattice suggests that the l-loop contribution to the effective action may 
be expanded in an asymptotic series of the form 



Close to the continuum limit, all terms with m > can be neglected. We more- 
over expect that the logarithmically divergent terms cancel after the gauge 
coupling go is renormalized, i.e. after we substitute 



00 




(4.20) 



m=0 n=0 



50 = aL + zi{ai^)gtt + . . . 



22 N 



ln(a/^). 



(4.21) 



3 167r2 
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where g\at is a renormalized coupling and /j, some normalization mass. This 
is entirely analogous to what happened in the case of the dimcnsionally reg- 
ularized effective action. In particular, the resulting renormalized expansion 
in powers of gf^^ must coincide with eq.(3.37), if the gauge couplings giat and 
gus are properly matched. 

To one-loop order and for the background fields considered, all these state- 
ments will be confirmed explicitly, and so we are confident that the dependence 
of the effective action on the lattice spacing is indeed as described. 

.^.5 0(a) improved action 

Numerical simulations of the functional integral (4.5) are limited to lattice 

spacings a that are not too small compared to the scales set by the background 
field. The cutoff dependent contributions to the effective action (the terms of 
order a and higher) may hence have a non-negligible influence on the outcome 
of such calculations. In that instance it may be desirable to choose an improved 
lattice action so as to cancel the dominant cutoff effects. 

This idea has previously been worked out for scalar theories [50-52] and 
pure gauge theories on lattices without boundaries [53-58] . The general princi- 
ple is that the cutoff effects of order a(lna)", a^(lna)", etc. can be successively 
removed by adding a linear combination of local counterterms to the action, 
with increasing dimensions and properly adjusted coefficients. This may be 
viewed as an extension of the renormalization procedure to the level of irrel- 
evant operators. We thus expect that the counterterms needed to improve 
the Schrodinger functional come in two forms, the boundary and the volume 
terms (cf. sect. 2). 

At order a there are two possible counterterms. They are obtained by 
summing any local lattice expression for the fields 

aUv{FokFok} and aHi {FkiFu} (4.22) 

over the = and = T hyper-planes. A simple 0(a) improved lattice 
action is thus given by eq.(4.4), where 

^ ^Cs(gio) if p is a spatial plaquette at x'^ = or = T, 

Ct(S'o) if p is a time-like plaquette attached to a (4.23) 
boundary plane, 

and w{p) = 1 in all other cases. The coefficients Cs{go) and ct{go) multiplying 
the boundary plaquettes allow us to monitor the size of the counterterms. 



w{p) = < 
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They are to be adjusted so as to achieve the desired improvement of the 
theory. 

Ideally one would like to determine Cs{go) and Ct{go) through numerical 
simulations of the Schrodinger functional. While this is not an impossible task, 
it is surely demanding, since one needs precise data on a range of lattices to 
be able to isolate the cutoff effects properly. In perturbation theory we have 

c,{g,) = cf)+c^})gl + ..., (4.24) 

Ct{go)=cf^ +c^^^9l + .--, (4.25) 

where the coefficients c'P and cf ^ can be extracted from the ^-loop contribution 
to the effective action. In particular, at tree level we only need to work out the 
order a term in the small a expansion of ro[-B] for two independent choices of 
B, say a constant Abelian field and a self-dual field. As a result one obtains 

4°) = cf^ = 1 for ah AT. (4.26) 

On the basis of our calculations to one-loop order (sect. 7), we have moreover 
been able show that 

c(i) = -0.166(1) and c\^^ = -0.0543(5) for N = 2. (4.27) 

Note that in the SU(2) theory the crossover from strong to weak coupling 
behaviour occurs around g^ = 2 and simulations nowadays are performed at 
couplings close to 1.5. The one- loop correction (4.27) is thus reasonably small. 

We finally mention that the improved action coincides with the Wilson 
action to lowest order of g^. The latter is hence 0(a) improved at tree level, 
and this implies that all coefficients r/^„[S] in the expansion (4.20) with 
m = 1 and n = I vanish. 
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5. Abelian Background Fields 



As explained in sect. 2, a large variety of running couplings can be defined 
by differentiating the effective action T[B] with respect to some parameter of 
the background field B. A careful choice of B is however necessary, if one 
attempts to study the scaling properties of the coupling through numerical 
simulations of the lattice Schrodinger functional. In particular, since one can- 
not afford to make the lattice spacing arbitrarily small, a background field is 
required for which the lattice corrections to the effective action are tolerable. 
It is our experience that constant Abelian fields are optimal in this respect 
and so we discuss their properties here. 

In the following we allow the space-time dimensionality D to assume any 
integer value greater or equal to 2. For the lattice action we take Wilson's 
action (subsect. 4.2) multiplied with a^""^. 

5.1 Definition 

The boundary values of the background fields considered in this section are 
spatially constant and diagonal. In other words, we have 



Cfc 



L 



fei 



0, 





'k2 



\ 



\ 




(5.1) 



where the angles (pf^ satisfy 



N 



'ka 



for all k. 



(5.2) 



The other boundary field C" is defined similarly, with (f)/^^ replaced by a second 
set of angles 4>'f^^. 

Gauge fields of this type have previously occurred in the literature and 
are referred to as "torons" (sec e.g. refs. [59,60]). Locally they are pure gauge 
configurations, i.e. the only gauge invariant information is contained in the 
parallel transporters for closed curves winding "around the world". In the 

field C these are products of the matrices exp{LCfc} and their inverses. In 
particular, they only depend on the phase factors e*'^'=« and the angular char- 
acter of the parameters (p/.^ thus becomes apparent. 



30 



An obvious solution of the field equations with the boundary values spec- 
ified above is 



Bo = 0, Bk= [x'^C'k + (T - /T. 
For the associated field tensor one obtains 

Gok = [C'k — C^.] /T, Gki = 0, 



(5.3) 



(5.4) 



and the action is given by 



S[B] 



L 



D-3 



D-1 N 

Yl Yl ("^fec 



^0 fc=l a=l 



(5.5) 



This solution — a constant colour electric field — is not the only solution of 
the Yang-Mills field equations with the required boundary values. It is in fact 
quite obvious that other solutions must exist, because the action (5.5) does 
not reflect the expected periodicity in the angles The stability of the field 
(5.3) thus needs to be discussed and we shall come back to this problem in 
subsect. 5.2. 

The hnk field 

y(.T, ;u) = exp {ai?^(a;)} (5-6) 

(with B as above) is a candidate for the induced background field on the 
lattice. V has the required boundary values and one may easily verify that it 
is a solution of the lattice field equations (4.12). Note that the action 



TL 



D-l 



9l 



D-l N 

EE 

k=l a=l 



sm 



2TL 



{<P'ka - <Pka) 



(5.7) 



coincides with the continuum action (5.5) up to terms of order (rather than 
a?'). Constant Abelian fields thus show their quality of leading to minimal 
cutoff effects already at the classical level. 
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5.2 Stability theorem 

Before the lattice field (5.6) can be accepted as the induced background field 
for the chosen boundary vahics, we must show that it is a configuration with 
least action and that any other field with the same action is gauge equivalent 
to V. Such a proof can in fact be given if (pf,^ and (p'f.^, are restricted to a 
certain bounded domain. 

To specify this region it is useful to arrange (?!)^^ (and similarly (^^^) into 
D — 1 colour vectors <I>^ in the obvious way. We say that a vector $ = 
(<^i, . . . , (f)j^) is in the fundamental domain if 

N 

(t>^<(t)2< ...<(t>N, I'/'a-'^/sl <27r, Y,^a = 0- (5-8) 

a=l 

This is a bounded convex set, which has a certain group theoretical signifi- 
cance. Namely, if u G S\J{N) has pairwise different eigenvalues A„, there exists 
a unique ordering of the eigenvalues and a unique vector $ in the fundamental 
domain, such that = e*"^" . 

Our result on the stability of the background fields considered in this 
section is now summarized by 

Theorem 1. Suppose the angle vectors and are in the fundamental 
domain for all k = I, . . . , D — 1. Let V be the associated background field 
(5.6) on a lattice with 

TL/a^ >{N- l)7r^ max{l, A^/16}. (5.9) 

The action of any other lattice gauge held U with the same boundary values 
then satisfies S[U] > S[V], where the equality holds if and only ifU is gauge 
equivalent to V. 

The proof of the theorem is quite involved and is deferred to appendix B. The 
condition (5.9) is of a technical nature and is perhaps not really needed for 
the validity of the theorem. It is, in any case, insignificant for < 3 and the 
lattices of interest. 

There is, of course, a continuum version of the theorem, provided one 
is willing to restrict attention to the category of differentiable fields. The 
minimum property of the solution (5.3) is in fact a simple corollary of theorem 
1, since any differentiable gauge field can be approximated arbitrarily well by 
lattice fields. 
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5.3 Definition of a running coupling for the SU(2) theory 

We are now ready to introduce the running coupling g^ (L) which will be used 
in our finite size scaling study of the SU(2) theory in four dimensions and 
which will, therefore, be in the focus of interest for the rest of this paper. 

Following the general scheme described in subsect. 2.6, we choose the solu- 
tion (5.3) for the background field B and set T = L. B depends on altogether 
6 parameters, say (f)i^-^ and (^^j. A symmetrical one-parameter submanifold is 

(l^kl = -11, ^fcl = -TT + 77, (5.10) 

and the coupling is now given by eq.(2.43). 

We should of course make sure that the premises of theorem 1 are satisfied 
so that the stability of the chosen background field is guaranteed. This is the 
case if 

< ?7 < vr and L/a> 4. (5.11) 

In the numerical work we take r] = 7r/4, which is half-way between the zero 
action point r] = ti/2 and the boundary of the stability interval. All these 
details are quite arbitrary, but we emphasize (once more) that they are in no 
way of fundamental importance for the scaling analysis we have in mind. 



6. Fixing the Gauge in the Lattice Theory 

The problem of gauge fixing has already been addressed in sect. 3 on 
a somewhat formal level. In particular, a truly convincing derivation of the 
boundary conditions on the Faddeev-Popov ghosts and the time component 
of the gauge field could not be given and so we feel that it is worthwhile to go 
through the gauge fixing procedure once more, in the mathematically rigorous 
framework of the lattice theory. We closely follow the general scheme discussed 
in subsect. 3.5 of ref.[47], where further details and a proof of the basic lemma 
can be found. 

In this section the background field is assumed to be an Abelian field as 
described in sect. 5, with angle vectors and in the fundamental domain. 
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6. 1 Group of gauge transformations 

The lattice Schrodinger functional (4.5) is invariant under all gauge transfor- 
mations U which leave the boundary fields W and W intact. According 
to lemma 1 of appendix B, this condition is satisfied if and only if the bound- 
ary values of the gauge function 0,(x) at x*^ = and x*^ = T are constant and 
diagonal. The group of all these transformation functions is denoted by Q. 

The constant diagonal gauge functions Q,(x) form a subgroup of Q iso- 
morphic to the Cartan subgroup Cat of SU{N). Such transformations act 
trivially on the background field and there are actually no further transforma- 
tions with this property (cf. appendix B). The gauge directions in the space 
of infinitesimal deformations of the background field are thus generated by 

G = G/Cn, (6.1) 

and so it is this set of transformations which needs to be fixed. The group Cjv 
then survives as a global symmetry of the theory. 

In the following it is convenient to identify Q with the group of all trans- 
formations CI e Q that are equal to 1 at x° = T. Prom a purely mathematical 
point of view, we then have the canonical situation, with a compact gauge 
group acting freely on the space of fields integrated over. The discussion in 
ref.[47] thus carries over literally. 

6.2 Spaces of infinitesimal fields 

The Lie algebra JC of G consists of all fields uj{x) such that the infinitesimal 
transformation 

n{x) = 1 - gou{x) + 0(^2) (6.2) 

belongs to Q. In other words, u> must be a spatially periodic lattice field which 
takes values in the Lie algebra of SU(A/") and which satisfies the boundary 
conditions 

where k is constant and diagonal. 

A linear space Ti. of lattice vector fields g^^ {x) may be introduced similarly. 
We here require that 

U{x,ii) = {l + goaq^{x) + Oigl)}V{x,ii) (6.4) 
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is a valid infinitesimal deformation of the background field. In particular, to 

guarantee that U has the same boundary values as V , the spatial components 
qk{x) must be equal to zero at x'^ = and = T. 

We shall later find it useful to have a scalar product on these spaces at 
our disposal. The obvious choice for the product of two vector fields is 

{q, r) = -2a^ ^ tr {g^(x)r^(x)} , (6.5) 

and the product on £. is defined similarly. 

Another notational item are the covariant lattice derivatives and D^. 
These operators act on functions /(x) with values in the Lie algebra of SU(iV) 
and are given by 

D^fix) = i [V{x, fi)f{x + af,)V{x, - f{x)] , (6.6) 

D;f{x) = I [fix] - V{x - aft, ii)-^f{x - aii)V{x - afi, m)] . (6.7) 

They only make sense if / is defined at the lattice points referred to and should 
thus be used with care. 

6. 3 Gauge fixing function 

The gauge fixing term to be added to the Wilson action is the square of 
a suitable gauge fixing function F. As explained in ref.[47], F must be a 
regular mapping from the space of fields integrated over to the Lie algebra C. 
For perturbation theory it is actually sufficient to specify F on an arbitrarily 
small but finite neighborhood of the background field V . The fields U in such 
a neighborhood may be parametrized by 

U{x,n) = exjp {goaQtiix)} V{x,n), (6.8) 

where q e H and say \\q\\ < e. The simplest possibility then is to choose the 
gauge fixing function to be a linear mapping from TC to C 

To be able to write it in a compact form we introduce the operator 

d: C<-^n, {duj)^ (x) = D^uj{x). (6.9) 

The covariant derivative occurring here is well-defined, for all relevant lattice 
points X and directions ji, and it is also easy to verify that duj has the proper 
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boundary values. The significance of d becomes apparent if we note that the 
gauge directions in the space of infinitesimal deformations of the background 
field are precisely the modes of the form duj. As already mentioned, the group 
Q acts freely on the space of gauge fields, and so it is not surprising that the 
kernel of d turns out to be trivial, i.e. dw = implies a; = 0. 

The basic property the gauge fixing function should have is that it does not 
vanish on the gauge modes dm. It is quite obvious now that this requirement 
(and all the other conditions listed in ref.[47]) will be fulfilled if we choose 

F{U) = d*q, (6.10) 

where d* is minus the adjoint of d. That is, d* maps any vector field q & TC 
onto an element of C such that 

{d*q,u!) = —{q,du!) for all u E C (6-11) 
This property implies 

d*q{x) = D^q^ix) if < x° < T, (6.12) 
and for the boundary values at a;° = and a;° = T one obtains 



[d*q{ 



X 



a/3 



UaVL^)Ey[QoiO,y)Uf, ifa = /3andxO = 0, ^^^^^ 
1 otherwise. 

Note that these are indeed as required for a function contained in £. 
6.4 Gauge fixed functional integral 

Now that we have specified the gauge fixing function F, the gauge fixed form 
of the Schrodinger functional (4.5) is obtained almost mechanically, following 
the steps described in ref.[47]. There is no new clement involved here and so 
we merely introduce the necessary notations and state the final result. 

Our choice of gauge fixing function corresponds to the gauge fixing term 

S^f[B,q] = ^{d*q,d*q). (6.14) 

In view of eq.(6.12), this may be regarded as a lattice version of the background 
gauge. The associated Faddeev-Popov ghosts c and c are in all respects like 
infinitesimal gauge transformations except that they are fermi fields. That is. 
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if we choose a complete set of functions in C, the ghost fields are equal to 
the general linear combination of these basis elements, with anti-commuting 
coefficients that generate a Grassmann algebra. 
The action of the Faddeev-Popov fields is 

Spp [5, g,c,c] = -(c,(i*(5cg), (6.15) 

where 5cq denotes the first order variation of q under the gauge transformation 
generated by c. To second order in we have 

6cq^l = D^c + 5oAd c 

+ \gQakdq^ + ^{g^aAdq^f + . . . D^c (6.16) 

(no sum over is implied here). Note that 5cq is a vector field with the correct 
boundary values. The action of d* in eq.(6.15) is therefore well-defined. 
The gauge fixed form of the Schrodinger functional is now given by 

e-r[B] ^ j j^p^ j D[c]D[c]e-^'°*-'[-^'«''='^l, (6.17) 

'S'total 

[B, q, c, c] = S[U] + S^i[B, q] + Sfp[B, q, c, c], (6.18) 

where it is understood that U and q are related by eq.(6.8). The first integral 

in eq.(6.17) is restricted to a small neighborhood of the background field V, 
i.e. wc are neglecting terms that arc exponentially small in the coupling. The 
integration over the ghost fields is the usual one, resulting in the determinant 
of the Faddeev-Popov operator. 

As in the continuum theory, the gauge fixed functional integral is the 
starting point from which the perturbation expansion of the effective action is 
obtained. Noting 

D[U]=B[q]{l + 0{g^o)}, (6-19) 

Stot.i[B,q,c,c] = S[V] + i {q,Aiq) + (c, Aqc) + 0(5o), (6.20) 
we get, for the first two terms in eq.(2.38), 

To[B]=glS[V], (6.21) 
ri[5] = ilndet Ai - Indet Aq. (6.22) 
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Ao and Ai are symmetric operators acting on functions in L and H respec- 
tively. We have already used the symbols Aq and Ai for the corresponding 
operators in the continuum theory, but this will not lead any confusion in the 
following. 

Prom the above we infer that 

Ao = -d*d, (6.23) 
Ai = A'l - Aodd*, (6.24) 

where A'l is obtained by expanding the Wilson action S\U\ to second order 
in q. With the help of the ★ product notation introduced in appendix A, the 
result of that computation may be written in the compact form 

A'ig^(x) = J^jcosh (a^G^,) * \-BlB^q^{x) + BlD^q,{x)\ (6.25) 

- oT'^ sinh (a^G^^) ★ ^q^(x) + a (£>* + D^) q^{x) + a'^D*Dnqu{x)] }. 

It should be emphasized that this formula is only valid for constant Abelian 
background fields, the case of immediate interest to us. We have also assumed 
that the link {x,x + afi) is contained in the interior of the lattice, viz. 

< a;° < r if /Lt = and < x° < r if > 0. (6.26) 

All degrees of freedom of the fields in H reside on these links and the operator 
A'l is thus completely specified by eq.(6.25). 

6.5 Boundary conditions 

The lattice quantum field q and the ghost fields c and c do not satisfy any 
boundary conditions beyond those implicit in the definition of 7i and C. In 
particular, the time component qo{x) is defined at all points x with < x^ <T 
and is unconstrained. 

However, after passing to the continuum limit, the situation regarding 
boundary conditions may be quite different. In a free scalar theory, for ex- 
ample, the propagator on a lattice with free boundary conditions converges to 
a Green function which satisfies Neumann boundary conditions. We should 
thus be prepared to find that the cigcnfunctions of the lattice operators Aq 
and Ai satisfy additional boundary conditions in the limit a ^ 0. 
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To be absolutely clear, we are not suggesting that the lattice theory should 
be amended in some way. Our aim rather is to determine which boundary 
conditions must be chosen when one attempts to formulate the theory directly 
in the continuum, using dimensional regularization for example. By comparing 
with the lattice theory this question can now be decided, because all cutoff 
prescriptions should yield identical results at tree level. 

In the case of the lattice operator Ai it is not immediately obvious how 
to pass to the continuum limit, because the definition of the gauge term dd* is 
slightly non-uniform at the boundary [cf. eqs.(6.12) and (6.13)]. This difficulty 
can be removed by choosing a better notation. To this end we extend the time 
component qo{x) of the lattice field to all points with x° = —a and = T. 
Its values there are chosen such that 

DSqo{x) = d*q{x) at = and = T. (6.27) 

No additional degrees of freedom are thus introduced and the extension should 

be regarded purely matter of convention. 

The operator dd* is now given by the simple expression 

dd*q^{x) = D^D*q,{x), (6.28) 

for all X and /i in the range (6.26). Furthermore, we may interpret eq.(6.27) as 
a boundary condition on qg. At x^ = T, for example, the requirement simply 
is that Doqo{x) = 0. At the other end of the lattice qo is again required to 
satisfy Neumann boundary conditions with the exception of the spatially con- 
stant diagonal modes, which must vanish at x^ = —a and thus satisfy Dirichlet 
boundary conditions. This latter complication depends on our choice of back- 
ground field and is absent for the irreducible background fields considered in 
sect. 3. 

In the present formulation one may straightforwardly pass to the contin- 
uum limit. In particular, taking a ^ in cqs.(6.25) and (6.28) leads to the 
continuum expression (3.18) for Ai and the boundary conditions stated above 
carry over literally. To remove all doubts about this procedure, we have verified 
in a number of cases, by exact numerical computation, that the eigenfunctions 
of the lattice operator converge to smooth functions in the continuum limit 
and that these indeed satisfy the expected boundary conditions. 

We thus conclude that the proper boundary conditions on the quantum 
field q in the continuum theory are a mixture of Dirichlet and Neumann bound- 



39 



ary conditions as described above. A similar analysis applies to the Faddeev- 
Popov fields. The result here is that besides the boundary conditions already 
present on the lattice, the spatially constant diagonal modes satisfy Neumann 
boundary conditions at = 0. This is again special to our choice of back- 
ground field in this section, i.e. for irreducible boundary values C and C one 
would end up with Dirichlet boundary conditions for all modes. 



7. First Order Correction to the Running Coupling [SU(2) Theory] 

We now set N = 2 and consider the one-parameter family of background 
fields defined in subsect. 5.3. Our aim is to compute the associated running 
coupling g'^ (L) to one-loop order of perturbation theory. The calculation will 

be performed in the framework of the lattice theory, but as a check we have 
also worked out the relevant determinants in the continuum theory, using 
dimensional regularization. Some of the details involved in that calculation 
will be sketched at the end of this section. 

7.1 Preliminaries 

In the following we choose lattice units and thus set a = 1. From the pertur- 
bation expansion of the effective action we deduce that 

f{L) =gl+ m^{L)gt + m2{L)gl + ... (7.1) 

The one-loop coefficient is a sum of two contributions, 

mi = ho - \hi, (7.2) 

one from the ghost and the other from the vector fields. Explicitly, they are 
given by 

/t, = — IndetA, /— ro[5], s = 0,l. (7.3) 
By differentiating the classical action (5.7) one obtains 

1 



^To[B] = -2AL' sin 



(7.4) 
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and we are thus left with the problem to compute the determinants of Aq 

and Ai for the background fields of interest. We do not expect that this can 
be done analytically, but as will be shown in the following subsections, it is 
possible to evaluate the determinants numerically, for a useful range of L and 
to essentially any desired precision. 

7.2 Symmetries 

Aq is a symmetric operator which acts in a real vector space. In particular, 
there exists a complete set of eigenfunctions and the determinant of Aq is the 
product of all eigenvalues. 

In the following we consider Aq to be a linear operator in the complex 
space C ® iC. The eigenfunctions and eigenvalues arc not affected by this 
extension and so it is clear that also the determinant is the same as before. 
The reason for going to the complexified space is that now we can more easily 
pass to a basis which is adapted to the symmetries of the problem. 

To exploit the invariance of Aq under constant diagonal gauge transfor- 
mations, we choose the standard basis of SU(2) generators [eq.(A.4)] and 
define = it iT^ . We then consider the decomposition 



where Cq and C± are the subspaces of all fields that are proportional to 
and respectively. Constant diagonal gauge transformations reduce to a 
multiplication by a phase factor on and so we expect that these spaces are 
invariant under the action of Aq. 

To show this explicitly, we introduce the corresponding subspaces TIq and 
7i± of complex vector fields and note that the operator d maps any function 
a; € onto an element of Tier- More precisely, we have 



It is quite obvious now that d* maps H.^ into and so we conclude that Aq 
leaves invariant. 



(7.5) 




(7.6) 



where the phase (5 is given by 



^(t) = -(2/L2)[ryL + (7r-2r?)t]. 



(7.7) 
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The operator Ai is similarly reduced by the subspaces Tier- This is an 
immediate consequence of the identities 

cosh Gofc * = cos(7/2)T", 7 = {2/L'){'k - 2r?), (7.8) 
sinhGofe*^^ = 0, (7.9) 
sinhGofe*r=^ = ^isin(7/2)r=^, (7.10) 

and the discussion above of the covariant derivatives. In particular, we have 

det Ai = det Ai|^^ det Ai|^_ det Ai|^^ , (7.11) 

and the analogous factorization holds for the determinant of Aq. Note that 
the last two factors in eq.(7.11) are equal. 

A further factorization of the determinants is obtained if wc take into 
account that Aq and Ai are invariant under spatial translations. The eigen- 
functions of this group of symmetries are the plane waves e*P^ where 

p = 27rn/L, nu e -L/2 < Uk < L/2. (7.12) 

We are thus led to introduce the subspaces Caip) and Tlaip) functions 
which are proportional to e*^^ and have no other dependence on x. Since the 
momentum p is conserved under the action of A 1, it follows that 

detAi|^^=ndetAi|^^(p), (7.13) 
p 

and the same factorization holds in the case of the operator Aq. 

The simplification which has thus been achieved is substantial. In each 
symmetry sector the spatial dependence of the wave functions and also the 
SU(2) degrees of freedom are completely fixed. The problem is, therefore, 
effectively reduced to computing the determinants of a set of finite difference 
operators in one dimension. 
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1.3 Computation of ho 

Prom the discussion above it is evident that Aq is independent of the back- 
ground field on the subspace Cq. The contribution from this sector hence 
drops out and we are left with 

- Indet Ao = 2^ - Indet Ao|^^(p) . (7.14) 

So let us consider the action of Aq in the subspace C+{p) in some more detail. 
The functions in this symmetry sector are of the form 

a;(x) =V'(a;°)e'P^T+, (7.15) 

where ■0(i), < t < T, is a complex function with 

V'(O) = ij{T) = (7.16) 

[cf. eq.(6.3); the constant k vanishes because T"*" is off-diagonal]. 

On this set of functions Aq reduces to an ordinary second order difference 
operator, 

AoV'(t) = Ait)i;it + 1) + B{t)i;it) + C(i)V'(i - 1), (7.17) 
with coefficients given by 

A{t) = Cit) = -1, (7.18) 

3 

B(t) = 8-2^cos[pfe+/3(t)]. (7.19) 

k=l 

The determinant of such operators can be computed by solving a simple re- 
cursion relation. This is a known trick which is discussed in some generality 
in appendix C. In the present case we need to solve 

Ao^(t) =0, < t < T, (7.20) 

starting from the initial values 



V^(o) = 0, = 1- 
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(7.21) 



The determinant is then given by 



(7.22) 



Note that the calculational effort to solve eq.(7.20) is proportional to T, while 
one usually needs of order arithmetic operations to evaluate the determi- 
nant of a T X r matrix. 

To compute the derivative of the determinant with respect to ij, the best 
way to proceed is to derive a recursion for dip{t)/dr] by differentiating eq.(7.20). 
Both recursion relations are trivial to program, and so one is able to compute 
the coefficient for say L < 32 with a negligible amount of computer time. 

7.4 Computation of hi 

For notational convenience we always assume that the field qo{x) is extended 
one step beyond the lattice to all points with x° = — 1 and x° = T, as dis- 
cussed in subsect. 6.5. The boundary conditions on then are a mixture of 
Dirichlet and Neumann boundary conditions and the action of Ai is given by 
eqs.(6.24),(6.25) and (6.28). 

We first consider the contribution of the sector 7^o(p)- The fields in this 
subspace may be written as 



where tp^it) is some complex vector field. The extra phase factors in these 
equations have no particular significance except that they lead to simpler 
expressions later on. 

The action of the operator Ai on wave functions of the above type is of 
the general form 



AiV'^(t) = A^.{t)Mt + 1) + B^,,{t)Mt) + C^At)Mt - 1)- (7-25) 



It is not difficult to work out the coefficient matrices A, B and C explicitly 
[cf. appendix D]. 

For all momenta p 7^ the boundary conditions on tpnii) 



qk{x) 



V'o(x°)(-z)e'P"r3, 
V^fc(x°)eiP'=e*P'^^^ 



(7.23) 



(7.24) 



d*ipo{t) = i/jkit) = at t = and i = T. 



(7.26) 
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d* here denotes the backward difference, defined in the usual way. If p = 
it is V'o(— 1) instead of 5*V'o(0) which is required to vanish, while all other 
boundary conditions are as above. 

The determinant of Ai in this latter sector can actually be computed 
immediately. The matrices A, B and C assume a particularly simple form in 
this case, leading to a decoupling of the -00 and ■0^ components. Furthermore, 
only the spatial modes make a field dependent contribution and one quickly 
finds that 

1^ Indet Ail^^(o) = 3(T - Incos |. (7.27) 

This is nice because now we are left with the p 7^ subspaces, where the 
boundary conditions are simply given by eq.(7.26). 

To compute the associated determinants, we proceed essentially as in the 
case of the Faddeev-Popov operator. We first construct all solutions of 

AiV'o(t) = 0, < i < T, (7.28) 
Ai^fe(t) = 0, < t < T, (7.29) 

with 3*V'o(0) = V'fc(O) = 0. For specified initial values, V'o(O) and dtpkiO), 
there exists exactly one solution. It is obtained recursively, first by solving 
eq.(7.28) at t = 0, then eq.(7.29) at t = 1, and so on. 

It is useful to group the initial data in a vector with components 

vo = MO), Vk = dMO)- (7.30) 
After completing the recursion we may compute the final values 

wo = d*MT), Wk = MT). (7.31) 

It is evident that these are related linearly to the initial values, i.e. there exists 
a matrix M such that 

Wf, = Mf,^v„. (7.32) 
The determinant of Ai is then given by 

det Ail^^(p) = (cos 1)'^"' detM, (7.33) 

as one may show by adapting the arguments of appendix C to the case at 
hand. 
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Table 1. One-loop coefiicient mi{L) and remainder e{L) at r/ = 7r/4 



L 


mi(L) 


siL) 


L 


mi(L) 


SiL) 


6 


0.35422030 


0.00352 


20 


0.47541031 


0.00024 


8 


0.38370429 


0.00177 


22 


0.48471069 


0.00020 


10 


0.40644167 


0.00107 


24 


0.49316916 


0.00017 


12 


0.42483107 


0.00072 


26 


0.50092488 


0.00014 


14 


0.44023578 


0.00052 


28 


0.50808526 


0.00012 


16 


0.45347745 


0.00039 


30 


0.51473490 


0.00010 


18 


0.46508332 


0.00031 


32 


0.52094162 


0.00009 



Let us now turn to the subspaces H+ (p) . The wave functions here are of 
the form 

qo{x)=M^'){-iy'"'T+, (7.34) 
= Vfc(^°)e*t^''+'^^''"^]e^P"T+, (7.35) 

and the action of Ai is again given by eq.(7.25), with coefficients A, B and 
C listed in appendix D. Furthermore, the boundary conditions are as before 
[eq.(7.26)], and so it is clear that the determinant of Ai in this sector can be 
computed following the lines discussed above. In particular, if M is the matrix 
connecting initial and final values, we have 

det Ail^^(p) =A^detM, (7.36) 

for all momenta p. 
7.5 Results 

We now set = 7r/4 and compute the one- loop correction to the running 
coupling using the techniques described above. Some of our results are listed 
in table 1. The calculations were done so as to guarantee a final precision 
of at least 20 decimal places. In particular, all digits quoted in table 1 are 
significant. 

In lattice units the continuum limit is approached by taking L to infinity. 
As discussed in subsects. 4.4 and 4.5 we expect that an asymptotic expansion 
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of the form 

oo 

mi(L) ~ V(r„ + s„lnL)/L" (7.37) 

L^oo '-^ 
n=0 

holds, with 

so = Y^ and si=0. (7.38) 

The first few coefficients r„ and s„ may be extracted from the data hsted 
in table 1 by any suitable extrapolation procedure. In particular, using the 
blocking method of ref. [58] , we were able to confirm, to four significant decimal 
places, that sq is given by eq.(7.38). This again shows that all ultra-violet di- 
vergences in the Schrodinger functional cancel after renormalizing the coupling 
in the usual way. 

If we now assume that sq is precisely given by eq.(7.38), the extrapolation 
of our data yields 



ro = 0.202349(3), (7.39) 
n = -0.1084(11). (7.40) 

As shown by table 1, the remainder 

e{L) = mi{L) - solnL - ro -ri/L, (7.41) 

with coefficients as given above, is rapidly decreasing over the range of L cov- 
ered. The higher order cutoff effects are thus comfortably small for our choice 
of background field. It must be emphasized, however, that other background 
fields, such as the constant self-dual field considered in sect. 2, do not fare as 
well in this respect. Partly this is because the strength of the Abelian field 
is constant while in the other cases the lattice must be fine enough to resolve 
the often appreciable time dependence of the field. 

After eliminating go in favour of the renormalized coupling (7iat introduced 
in subsect. 4.4, the continuum limit can be taken and one ends up with the 
expansion 

9L + --- (7.42) 

At this point we may pass to any other renormalization scheme defined in the 
continuum theory. In particular, if we choose the running coupling 5ms (-^) 



11 
12^ 



ln(/iL) -I- ro 
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the MS scheme of dimensional regularization as our new expansion parameter 
and use the known relation between g-iat and ^ms [62,9,63], we obtain 



(7.43) 



where the coefficient ci is given by 



11 



[In47r-7E - 1.61638(8)] . 



(7.44) 



ci = 



247r2 



For the corresponding A parameter ratio one finds 



A/Ams = 1.18378(5), 



(7.45) 



which shows that our coupling is very nearly equal to 5ms (-^) pertur- 
bative domain. 

The term proportional to ri in the expansion (7.37) is an order a cutoff 
effect. As discussed in subsect. 4.5 it can be removed by adding a combination 

of two boundary counterterms to the Wilson action, with coefficients Cg and 
Ct- Our result above, eq.(7.40), allows us to compute q to one-loop order 



To also determine Cg, the effective action of a further background field, 
with non-vanishing spatial field strength, must be computed. The constant 
self-dual fields are a convenient choice for this calculation, because the deter- 
minants of Aq and Ai again factorize after going to momentum space. With 
some additional work, we have thus been able to obtain the result quoted in 
subsect. 4.5. 

7.6 Computation of g^{L) using dimensional regularization 

It is possible to derive the expansion (7.43) directly in the continuum the- 
ory, starting from the dimensionally regularized Schrodinger functional. This 
provides an important check on our calculations on the lattice and gives us 
additional confidence that the Schrodinger functional is indeed a universal 
amplitude. 

In sect. 3 we have already obtained the one-loop effective action for a gen- 
eral irreducible background field, using dimensional regularization. The final 
formula, eq.(3.37), is also valid for Abelian background fields, if Aq and Ai 
are considered to be operators in the proper function spaces [cf. subsect. 6.5]. 
The non-trivial remaining problem then is to compute i^'(0|Ao) and (^'(0|Ai) 
or rather their derivatives with respect to r]. 



[cf. eq.(4.27)]. 



48 



In the following we adopt the notation of sect. 3. In particular, the symbol 
A denotes any one of the operators Aq or Ai. Prom the definition of the zeta 
function and the Seeley-DeWitt expansion we infer that 



C'(0|A) = 7E«o(A)+ 




4 



ao{A)e{l-t) 



(7.46) 



This integral is absolutely convergent and the lower end of the integration 
range may, therefore, be replaced by some positive cutoff 6 which is subse- 
quently sent to zero. It is then straightforward to show that 



-C'(0|A) = lim (7E + ln^)^ao(A) - E ' (7-47) 



where En, n = 0, 1, 2, . . ., are the eigenvalues of A. 

Our computational strategy now is as follows. We first calculate all eigen- 
values and their derivatives dEn/drj up to a certain level. For values of 
S that are not too small, the bracket in eq.(7.47) is then accurately given by 
restricting the spectral sum to these eigenvalues. Finally, the result is ex- 
trapolated to the limit 5 = 0, taking into account that the bracket has an 
asymptotic expansion in powers of 6^^^. 

The success of the method of course depends on our ability to compute 
sufficiently many eigenvalues of A. This goal is easy to reach. We first fix 
the symmetry properties of the wave functions as before, by choosing an ap- 
propriate basis of group generators and by going to momentum space. The 
range of momenta p is here unbounded, but since the eigenvalues are growing 
essentially like p^, only a finite subset of momenta needs to be considered. 

In every symmetry sector the problem is then reduced to finding the 
eigenvalues of a certain second order ordinary differential operator acting on 
functions with Dirichlet or Neumann boundary conditions. The reduced op- 
erators are actually so simple that the eigenfunctions can be given explicitly 
in terms of hypergeometric functions with a parameter to be adjusted. Al- 
ternatively, one may use exact numerical methods (Runge-Kutta integration 
and root finding algorithms) or one may set up a variational calculation with 
a suitable basis of plane wave states. 



CXD 



n=0 
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Proceeding either way we have been able, with a small amount of com- 
puter time, to calculate the bracket in eq.(7.47) for values of 6 in the range 
between say 0.01 x and 0.001 x L^. From there an extrapolation to 5 = 
is safely possible and our result for the coefficient ci was found to be in agree- 
ment with eq.(7.44) (the estimated precision after the extrapolation was about 
4 significant digits). 



8. Concluding Remarks 

The renormalizability of the Schrodinger functional opens the door to a 
new generation of scaling studies of pure gauge theories. We have in this paper 
prepared the ground for such an investigation. Numerical simulations of the 
SU(2) theory are now underway and first results are expected soon. 

It would be quite trivial to extend our calculations in sect. 7 of the running 
coupling to the physically more interesting case of the SU(3) theory. Before 
that one however needs to fix the parameters of the background field. In 
particular, some experience must be accumulated to determine which fields 
lead to an acceptably large signal in numerical simulations, while keeping the 
cutoff effects low. 

We also believe that an extension of our computations to the two-loop 
level is feasible, although this is clearly going to be hard work. 



Appendix A 

Our notational conventions are as follows. Lorentz indices /i, i/, . . . nor- 
mally run from to 3. In the context of dimensional rcgularization they run up 
to D — 1, the dimension of space. Since the metric is euclidean, it does not mat- 
ter in which position these indices appear. The totally anti-symmetric symbol 
^nvpa- is normalized such that eoi23 = 1- Latin indices k,l,... run from 1 to 3 
(or D — 1) and are used to label spatial vector components. Colour vectors in 
the fundamental representation of SU(A^) carry indices a, (3, . . . ranging from 
1 to N, while for vectors in the adjoint representation, latin indices a,b, . . . 
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running from 1 to — 1 are employed. Repeated indices are automatically 
summed over unless stated otherwise. 

The Lie algebra su(A^) of SU(A^) can be identified with the space of com- 
plex N X N matrices which satisfy 

= -X, tT{X} = 0, (A.l) 

where denotes the adjoint matrix of X and tr{X} = X^a is the trace of 
X. This is a real vector space of dimension A^^ — 1. A natural inner product 
is given by 

{X,Y) = -2tv{XY}, (A.2) 
and we may thus choose a basis T"-, a = 1,2,..., N'^ — 1, such that 

For N = 2, the standard basis is 

r« = ^, a = 1,2, 3, (A.4) 

where r" denote the Pauli matrices. With these conventions, the structure 
constants /"^^, defined through 

[r",r''] = /"^^T^ (A.5) 

are real and totally anti-symmetric under permutations of the indices. 

An SU(A^) gauge potential is a vector field A^{x) with values in the Lie 
algebra su(A^). With respect to a basis T'^ of group generators, we have 

A^{x) = A;ix)T\ (A.6) 

where the component fields v4^(x) are real. The components of other fields 
such as the gauge field tensor or the Faddeev-Popov ghost field are defined in 
the same way. 

The representation space of the adjoint representation of su(Ar) is the Lie 
algebra itself, i.e. the elements X of su{N) are represented by linear transfor- 
mations 

AdX: su(iV) ^ su(iV). (A.7) 
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Explicitly, AdX is defined through 

AdX{Y) = [X, Y] for all Y G su(iV). (A.8) 

With respect to a basis T", the associated matrix (AdX)"^ representing the 
transformation is given by 

AdX(r'') =^"(AdX)"^ (A.9) 

which is equivalent to 

(AdX)"^ = -/"^^X^ {AdX{Y)y = f''^''X''Y^, (A.IO) 

in terms of the structure constants. 

We finally define the ★ product of two N x N matrices M and X through 

M-kX = ^ {MX + XM"^) - ^ tr {MX + XM"^) . (A.ll) 

M + X is contained in su(A'^), for all X G su(A?^) and arbitrary M. 



Appendix B 

The proof of theorem 1 is elementary from the mathematical point of 
view, but it is lengthy and is therefore broken up into several steps. Our first 
goal will be to show that the theorem holds in all dimensions if it is true in 
two dimensions. After that we shall work out the consequences of the lattice 
field equations in two dimensions. We shall find that there is a discrete set 
of solutions and we are then left to prove that among these the configuration 
with least action is the background field V. 

Without further notice, we assume in the following that the angle vectors 
$^ and are in the fundamental domain (5.8) and that the bound (5.9) 
holds. The discussion will be solely concerned with the lattice theory and 
the boundary values of the gauge fields considered are always as specified in 
subsect. 5.1. Furthermore, the letter V is strictly reserved for the background 
field defined through eqs.(5.3) and (5.6). 
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B.l Remarks on gauge transformations 

The boundary fields W and W have an important stability property under 
gauge transformations which is expressed by 

Lemma 1. If A is a spatial gauge transformation which leaves W or W 
invariant, it must be equal to a constant diagonal matrix. 

Proof: Let us assume that it is W which is invariant under the action of A, 
viz. 

A(x)VF(x, A;)A(x + ak)"^ = ■^^(x, k). (B.l) 

The parallel transporter along any lattice curve at = is an ordered prod- 
ucts of the link variables and so must be invariant too. In particular, we may 

consider a closed curve starting at some point x and winding once around the 
world in the fc-direction. Since A is periodic, we then deduce that 

A(x) exp {LCk} A(x)-i = exp {LCk} . (B.2) 

The parallel transporter exp {LCk} is a diagonal matrix with eigenvalues = 
Q^<t>ka . These are all distinct (the angle vector is in the fundamental domain) 
and eq.(B.2) thus implies that A(x) must be diagonal. That it must also be 
constant, is now an immediate consequence of eq.(B.l). D 

A simple corollary of this lemma is that the background field V is invariant 
under a space-time gauge transformation Q, if and only if Q, is constant and 
diagonal. 

For any given lattice gauge field U we define 

@{x) = U{y,0)U{y + a6,0)U{y+2a6,0)...U{x-a6,0), y = (0,x), (B.3) 

which is just the parallel transporter along the time axis from x to the "base 
point" y at the boundary. @ is periodic in the space directions and may be 
regarded as a gauge transformation function. 

Lemma 2. IfU is known to be gauge equivalent to V, wc have V = . 
Proof: Let us define U = ?7®. From the definition (B.3) of 6 we infer that 

U{x, 0) = V{x, 0) = 1 for all x, (B.4) 
U{x, k) = V{x, k) = Ty(x, k) for all k and all x with x^ = 0. (B.5) 
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Furthermore, since U and V are gauge equivalent, there exists a gauge trans- 
formation Q. such that U = . 

Eq.(B.4) now imphes that 17 is time independent, and if we combine 
eq.(B.5) with lemma 1, one concludes that O must be constant and diagonal. 
It follows that = V and so we have proved the lemma. D 

B.2 Reduction to two dimensions 

A significant simplification is now achieved by 

Lemma 3. Let us assume that theorem 1 holds in two dimensions. Then it 
is true for all D >2. 

Proof: The action of any field U in D dimensions may be written in the form 

S[U] = Se[U] + Sm[U], (B.6) 

where Se and Sm arc the contributions of the time-like and spatial plaquettes, 
respectively. The magnetic part Sm is always greater or equal to zero and so 
we conclude that 

S[U] > Se[U]. (B.7) 
The electric part Se may be split up according to 

SE[U]=a''-^J2^l^U (B-8) 

V 

The sum here extends over all two-dimensional time-like planes V on the lattice 
and U\^ is the two-dimensional gauge field that one obtains by restricting U 
to V. 

We now make use of the premise of the lemma to show that 

Se[U] > ^ S[V\-p] = Se[V] = S[V]. (B.9) 

V 

When combined with cq.(B.7), this proves that S[U] > S[V] for all fields U. 

We still need to show that U is gauge equivalent to V if S[U] = S[V]. So 
let us assume that U is a minimal action configuration. Prom the above we 
then infer that the two-dimensional fields U\.p and V\.p are gauge equivalent, 
for all planes V. This does not immediately imply that U is gauge equivalent 
to V, because the gauge transformations which map U\-p onto V\-p need not 
coincide at the intersections of different planes. 
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At this point lemma 2 comes to rescue. It is vahd in any dimension and 
may be apphed to U\j, and V\-p (in place of U and V). We then conclude 
that 

[U\rf = V\v^ (B.IO) 

where Q is the transformation defined by eq.(B.3). Since Q{x) is independent 
of the plane V passing through x, we have thus shown that J7® = V. D 

B.3 Group theoretical lemmas 

To be able to work out the consequences of the lattice field equations later on, 
some technical preparation is needed. 

The eigenvalues Aq of any matrix u G SU(A'^) lie on the unit circle and 
their product is equal to 1. We say that u belongs to the set <S, if 

N 

A„=e^^-, where |xa| < 7r/2 and ^ Xa = 0. (B.ll) 

a=l 

It is obvious that S is an open neighborhood of the group identity. A useful 
criterion for a matrix m to be an element of S is 

Lemma 4. Any matrix u G S\J{N) with Retr(l — u) < min{l, 16/A?^} is 
contained in S. 

Proof: The eigenvalues Aq, of u may always be parametrized by a set of angles 
Xa such that 

N 

A„=e^^«, -7r<Xa<7r, ^Xa = 27rn, (B.12) 

a=l 

where n is an integer. Prom the bound 

N 

Retr(l -u) = ^(1 -cosxa) < 1 (B.13) 

a=l 

we conclude that |xa| < 7r/2. To prove that the sum of the angles Xa vanishes, 
we first note that 



l-cosx> (2x/7r)' if |xl<vr/2. 
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(B.14) 



The sequence of inequalities 



TV 



a=l 



<iV5^Xa < ^Retr(l-n) <47r2 (B.15) 

a=l 



is then easily established. It follows from these that |n| < 1, and since n must 
be an integer, we are left with n = as the only possible value. \3 
The significance of the set S is elucidated by 

Lemma 5. Any two matrices u,v & S with the same traceless anti-hermitean 
parts, viz. 

u — — -^tr (u — u^) = V — — -^tr (v — v^) , (B.16) 

are equal. 

Proof: Let A„ = c'''^" be the eigenvalues of u as in cq.(B.ll). The associated 
eigenvectors are also eigenvectors of u — u'^ , with eigenvalues 2i sin Xa ■ Since 
Xa 7^ Xl3 implies sinxo, ^ sin Xi3, the eigenspaces of u and u — must in 
fact coincide. The same statement applies to the other pair of matrices, v and 
V — . So if we take eq.(B.16) into account, it follows that all four matrices 
can be simultaneously diagonalized. 

Let us now choose a basis of simultaneous eigenvectors. The eigenvalues 
of V may be parametrized in the same way as those of u through a set of angles 
ipa- Eq.(B.16) then reduces to 

sin Xa - sin V'a = (B.17) 

where t is independent of a. In particular, if t > one infers that Xa > ipa for 
all a (recall that all angles are between — 7r/2 and 7r/2). This is impossible, 
however, because both sets of angles must sum up to zero. The same argument 
excludes t <0 and so we are left with Xa = i'a as the only acceptable solution 
of eq.(B.17). The matrices u and v are hence equal. D 

B.4 Proof of theorem 1 for D = 2 

Prom now on the discussion is restricted to the two-dimensional theory. To 
simplify the notation we set ^i^, = (f)^ and (t)'i^ = (f)'^ . Since the set of all lattice 
gauge fields is a compact manifold, the infimum of the action is attained by at 
least one configuration U. In the following we assume that U is such a field. 
We then need to show that U is gauge equivalent to V. 
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Lemma 6. For all x the plaquette matrix 

P{x) = U{x, 0)U{x + aO, l)U{x + al, 0)-'^U{x, 1)"^ (B.18) 
is contained in the set S. 

Proof: Since U is a minimal action configuration, we have 

S[U] < S[V] < {glTL)-' J2 - 4>S . (B.19) 

If we remember that $ and are in the fundamental domain, it is easy to 
prove that the angle differences tpa = (f^'a ~ 0a satisfy 

N 

|V'a-V'/3| <27r, ^Va = 0. (B.20) 

a=l 

Using these properties, the estimate 

N I ^ 

E ^« = ^ E (V'a - < 2(iV - l)7r2 (B.21) 

a=l a,l3=l 

may be derived. 

The action S[U] is a sum of non-negative contributions, one from each 
(unoriented) plaquette. Any one of these must be smaller than the right hand 
side in eq.(B.19) and thus, taking eq.(B.21) and the bound (5.9) into account, 
one deduces that 

Retr(l-P(x)) <min{l,16/iV} (B.22) 
for all X. Lemma 4 now tells us that P{x) is contained in S. \3 
Lemma 7. The plaquette field P{x) is covariantly constant, i.e. it satisfies 

P{x) = U{x, n)P{x + afi)U{x, n)-^ (B.23) 

for all X and directions 

Proof: Being a configuration with least action, U must be a solution of the 
lattice field equations (cf. subsect. 4.3). In two dimensions these are equivalent 
to the requirement that the unitary matrices 

u = P{x) and v = U{x,ii)P{x + a(i)U{x,ii)~'^ (B.24) 
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have the same traceless anti-hermitean parts, for all x and directions /x. Prom 
the above we know that u and v are contained in the set S and lemma 5 thus 
implies that u = v, as was to be shown. D 

The plaquette field P{x) may be regarded as a particular gauge transfor- 
mation function. Eq.(B.23) then simply says that U is invariant under this 
transformation. In particular, the boundary field W is left invariant and so, 
by lemma 1, we conclude that P{x) must be constant and diagonal at a;° = 0. 
In other words, we have 



(B.25) 



/e'^' ... \ 
e*>^2 ... 
^(^) 1x0=0= : : ■. : 

\ ... e'^"" J 
where the angles Xa be chosen such that 

N 

|Xa|<7r/2 and ^ Xa = (B.26) 

a=l 

(as a consequence of lemma 6) . 

Another important implication of eq.(B.23) is that the eigenvalues of P{x) 
are the same for all x. In particular, the action 

N 

ITT , , 

'^[^] = 7274E(l-^0«Xa) (B.27) 

is determined by the angles Xa- 

Lemma 8. There exists a permutation a and a set of integers Ua such that 

,2 



Xc 



In the special case where a is the identity and ria = for all a, the conhgura- 
tion U is gauge equivalent to V. 

Proof: We first pass to the temporal gauge by applying the gauge transforma- 
tion @ [eq.(B.3)]. The transformed field U = and the associated plaquette 
field P satisfy 

U{x, 0) = 1 for ah x, (B.29) 
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U{x, 1) = W{yi, 1), P{x) = P{x) at x° = 0. (B.30) 

Since P{x) is covariantly constant, we conclude that it must be equal to -P(O) 
for all X. 

Next we note that 

P{x) = U{x + aQ,l)U{x,l)-^, (B.31) 

which immediately leads to 

U{x, 1) = [P{Q)f''' W^(x, 1). (B.32) 

The matrices on the right hand side of this equation are diagonal, i.e. we have 
shown that U is an Abelian solution of the field equations similar to V . 

Let us now consider the parallel transporter for a loop winding once 
around the world at x^ = T. In the temporal gauge it may be computed 
directly using eq.(B.32), or we may note that it is just the gauge transform 
of the parallel transporter determined by the original field U. We thus obtain 
the relation 

e{x)cxp{LC[}e{x)-^ = [P(0)]^^/"%xp{LCi}. (B.33) 

The eigenvalues of the matrix on the left are the phase factors e"^« . Up to a 
possible reordering, they must be equal to the diagonal elements of the matrix 
on the other side of the equation and so we conclude that cq.(B.28) holds for 
some permutation a and some integers rioi- Furthermore, if a happens to be 
the identity and na = for all a, it is immediately clear from eq.(B.32) that 
U = V, i.e. U and V are gauge equivalent in this case. D 

The upshot then is that the possible minimal action configurations U are 
labelled by a permutation a and a set of integers 71^. Not all configurations 
{a, n) can occur, but only those for which the angles (B.28) satisfy the con- 
straints (B.26). These will be called admissible in the following. In particular, 
the trivial configuration, where a is the identity and = for all a, is 
admissible. Note that the set of admissible configurations {a, n) is finite. 

Taking eq.(B.27) and lemma 8 into account, the proof of theorem 1 is now 
completed by 
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Lemma 9. On the set of admissible configurations {cr,n), the minimum of 
the function 



N 



s{cr,n) = ^(l-cosx„) 



(B.34) 



a=l 



is not degenerate and the unique minimizing conGguration is the trivial one. 

Proof: Let us assume that (a, n) is a non-trivial admissible configuration. Wc 
then show that there exists another admissible configuration (a, n), such that 
s(a, n) > s(a, n). 

We first consider a configuration (cr, n) where not all n^s are equal to 
zero. Since {a, n) is admissible, the integers tIq, must add up to zero. It is, 
therefore, possible to find two indices a and /3 such that 



n„ > 1 and n^j < —1. 

Let us now define a new configuration (a, n) through 

' o-(^) if 7 = a, 
5(7) = } a{a) if 7 = /3, 
^ (7(7) otherwise. 



(B.35) 



(B.36) 



and 



' Ufj + 1 if 7 = a, 
Ua-l if 7 = /5> 
n-v otherwise. 



(B.37) 



The associated angles Xj ^^re equal to Xj with the exception of 



Xc 



Xl3 



TL 
TL L 



+ 27r (n^ + l)-(f>a 



^'aM + K - 1) - 



(B.38) 
(B.39) 



Note that a is just a transposition of a and /3 followed by a. 

It is not difficult to show that {^,n) is an admissible configuration. Fur- 
thermore, a little algebra yields 

s{a,n) - s{a,n) = (B.40) 
4 cos [i {xa + X/3 + Xa + Xp)] sin [i ixa - Xa)] siu iXa " X/?)] • 
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The arguments of all trigonometric functions in this expression are less than 
7r/2 in magnitude. In particular, the first factor is positive. Concerning the 
other two factors, we note that their arguments 



Xa Xc 



X« - X/3 



a 
TL L 

TL 



[27r + 



MP) 



+ 27r {jia - n/3 - 1) 



(B.41) 
(B.42) 



are positive, because the angle vectors $ and are in the fundamental domain 

and because the integers Ua and np satisfy the bounds (B.35). So we conclude 
that the action of the new configuration (cJ,n) is strictly lower than s(cj, n). 

Let us now consider a configuration (a, n) , where a is non-trivial but 
where all n^'s vanish. All configurations of this type are admissible. Since a 
is not the identity, there are two indices a and (3 such that 



a < P and cr{a) > a{P). 



(B.43) 



The new configuration {a,n) in the present case is given by eq.(B.36) and 
n-y = for all 7. Eq.(B.40) is then still valid and in view of 



Xct Xce 



Xa - X/3 



a 

TL 
TL 



MP) 



(B.44) 
(B.45) 



all factors on the right hand side are positive. It follows that s(cr, n) > s{a, n) 
and we have thus proved the lemma. D 



Appendix C 



In this appendix we derive a useful expression for the determinant of a 
general second order difference operator acting on a space of wave functions 
with Dirichlet boundary conditions. This result is needed in sect. 7 to compute 
the determinants of Ao and Ai in a constant Abelian background field. The 
basic idea is borrowed from Coleman's Erice lecture on the uses of instantons 
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where he treats the analogous case of a second order differential operator (see 
ref.[61], p.340). 

C.l Definitions 

The general second order difference operator A acts on complex wave functions 
V'(t) with n components, defined at integer values of t. Explicitly, A is given 

by 

AV'(t) = A{t)i^{t + 1) + B{t)^{t) + C{t)i^{t - 1), (C.l) 

where A, B and C are some complex n x n matrices depending on t. 

The wave functions which are defined for t = 0, 1,2, ... ,r and satisfy 
Dirichlet boundary conditions, 

^(t) = at t = and t = T, (C.2) 

form a vector space J- of dimension dy^ = n{T — 1). If ip{t) is an element of 
this space, eq.(C.l) is meaningful for < t < T and A may thus be regarded 
as a linear operator in J^. 

In the following we assume that A{t) is invertible for all t and that A is 
hermitean relative to some scalar product on J^. In particular, we take it for 
granted that there exists a complete set of eigenfunctions. 

C.2 Statement of result 

For any complex A the equation 

(A - A)^(t) =0, t>0, (C.3) 

has a unique solution with V'(O) = and prescribed initial value at i = 1. It 
can be computed recursively by first solving eq.(C.3) at i = 1, then at t = 2, 
and so on. After T — 1 steps one obtains ijj{T) which is in general not equal 
to zero. It is obvious, however, that ipiT) depends linearly on the initial value 
-0(1), i.e. there exists an n x n matrix M(A) such that 

^(T) = M(A)V(1). (C.4) 

This matrix is evidently determined through the coefficients A, B and C. It is 
comparatively easy to evaluate, requiring a computational effort proportional 
to n^T. 



62 



The formula alluded to above now reads 

T-1 

det (A - A) = det M(A) JJ det [-A{t)] , (C.5) 

t=i 

where A is here considered to be an operator in J^, as described above. In 
particular, the determinant of A is obtained if we set A = 0. 

C.3 Proof of eq. (C.5) 

It is easy to show that the matrix M(A) is a polynomial in A of degree T — 1. 
The leading term is 

M(A) = A^-^ {A{1)A{2) ...A{T- 1)}-' + 0(A^-2). (C.6) 

It follows from these remarks that 

T-1 

P(A) = det M(A) Yl det [-Ait)] (C.7) 
t=i 

is a polynomial in A too, with leading term equal to (— A)*^^. 

Let us now assume that /x is an eigenvalue of A with multiplicity k. If 
ijj e is one of the associated eigenfunctions, we have 

= ^(T) = Miii)i,{l). (C.8) 

The matrix M(|u) thus has a zero mode and we conclude that P{n) = 0. 
Actually, since there are k linearly independent eigenfunctions, we may choose 
a basis such that the first k columns of M{p) vanish. The multiplicity of the 
zero of -P(A) at A = is, therefore, greater or equal to k. 

The total number of eigenvalues of A, including multiplicities, is equal to 
djF- Since this is also the degree of P(A), it follows that this polynomial cannot 
have any further zeros besides the eigenvalues of A and that, moreover, the 
corresponding multiplicities must coincide. 

We have thus shown that -P(A) is the characteristic polynomial of the 
operator A, which is precisely the content of eq.(C.5). 
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Appendix D 



In this appendix the matrices Af^,^, B^^, and C^i, introduced in subsect. 7.4 
are given exphcitly. They are all real and satisfy 



B^^.it) = B,^{t), C^,{t) = A.^{t - 1), 



(D.l) 



as required for a symmetric operator. In the following we shall use the abbre- 
viations 



Pk = 2 sin [\pk] , 
Sk(t) = 2sm[\{pk + m)\ ■ 



(D.2) 
(D.3) 



We now list the independent elements of the matrices which describe the action 
of Ai in the HoCp) sector. 



-^00 (0 — ~-^0) 

Aki{t) = - cos - 5ku 



Ako{t) = 0, 



cos An 

2 



Pk, 



Booit) = 2Ao + cos -pjPj, 



Bki{t) 
Bok{t) 



T 1 
2 cos - +PjPj^ 5ki + (Ao - l)pkPi, 



cos An 

2 



Pk, 



(D.4) 
(D.5) 

(D.6) 
(D.7) 

(D.8) 

(D.9) 

(D.IO) 



The matrices describing the action of Ai in the H+(p) subspace are as follows. 

Aoo{t) = -Ao, (D.ll) 

Aki{t) = -Ski, (D.12) 

Aokit) = XoSkit + 1) - Skit), (D.13) 

Ako{t) = 0, (D.14) 
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